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RÉSULTATS EXPÉRIMENTAUX 
‘ET INTERPRÉTATION 


… Expériences sur le mercure, — Les différentes 
- expériences ont été faites sur des tubes électro- 

- niques remplis avec l’un des mélanges d’isotopes 
- suivants : 


I IL III 
(MERCURE NATUREL) 


ne. Isotopes pairs 69,8 % 20 0 MO 
199 417,0 % 73410 0, 0,8% 
201 13,200 


2500 CM88:2 70 


I. Isoropes PAIRS. — Étude de la forme de raie. 
— L'étude de la forme d’une raie de résonance 
magnétique a été étudiée en détail par Brossel et 
Bitter dans le cas du niveau 6 ?P, du mercure [1], 
[2], [3]. Nous allons généraliser cette étude au cas 
d’un niveau quelconque. : 

Supposons que, dans un régime stationnaire 
d’excitation, le nombre d’atomes dans le sous- 
niveau Zeeman m — p soit NV,, nous pouvons faci- 
lement calculer la grandeur Z, — 1,, grandeur qui 
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Résumé. — Dans un article précédent, nous avons décrit les expériences permettant de mettre 
en évidence la résonance magnétique des niveaux 6 %F,, 54° 6s? 6p et d’un niveau ?P, de l’atome 
de mercure excité par bombardement électronique. : 

Nous montrons que l’étude de la largeur des raies des isotopes pairs conduit à la durée de vie du 
niveau. Nous avons déterminé d’autre part les facteurs de Landé. Dans le cas des isotopes impairs 
Hg et ?2Hg, nous avons étudié l'effet Back-Goudsmit et déduit un ordre de grandeur des struc- 


Des résonances magnétiques ont été décelées par la même technique sur l’atome de sodium. 


Abstract. — In a previous paper experiments have been described in which magnetic resonances 
of levels 6 %F,, 5d° 6s? 6p D, and $P, of mercury were produced. 1 


Electron bombardment was 


It is shown that measurements of the width of the even isotope lines lead to the value of the 
Landé factors have also been determined. 

In the case of the odd isotopes Hg and Hg, the Bàack-Goudsmit effect of the 6 5F, level has 
been studied. An approximate value of the hfs intervals of this level has been obtained. 

The same technique has given positive results on some levels of the Na atom. 


nous est accessible grâce au montage en pont des 
photomultiplicateurs. On a : 


p=—m 


Tr — Is = D 


p=m 


No Cp (L) 


les €, sont les différents coefficients numériques 
calculables à partir des probabilités de transition. 

Envisageons une excitation par chocs électro- 
niques ne durant que le temps dt à l'instant 0, 
nm, atomes sont portés par seconde dans lPétat 
excité m — p. Au bout du temps {, par émission 
spontanée, le nombre d’atomes sera réduit à 
ny e * dt, puisque la durée de vie x est la même 
pour tous les sous-niveaux Zeeman. 

Si, dans ces conditions, nous établissons un 
champ de radiofréquence pendant le temps t, des 
atomes seront transférés des états m = p aux 
états m —q avec une probabilité P(F, p, q, 1) 
définie par aillleurs (*). La variation de popu- 
lation du niveau M — p sera : 

g= +m 


dnp = > 


= 


P(F, P; q; t)ng 


a=+m 
D PE opt om r 0e 
q=—m = 
(*) L’expression de Majorana, P(F, p, q,1), est donnée en 
particulier dans [2] et [3]. 
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soit, comme P(F, p, q,t) = P(F, q, p, t) : 
din = 


D PERS P; 4 !) (ra AT ny) etlr dt. 
q 


Si l'excitation électronique se produit pendant 
longtemps, la variation de population du niveau 
m = psera : 


Any = 2 1e P(F, p, q t) (na — np) eT #7 dé. 
ga V0 


D’après (1), le signal, lors de la résonance magné- 
tique, est : : 


s—= As — 1) =Y Any Cr 
D 


SOI 


D=+m = +m 


FE DEES 


D=— M q=—M 


CA 
LP p, @ 0 (na — mn) Cre#r de 
0 


la vitesse d’excitation se déduit du nombre AN, 
d’atomes dans un régime stationnaire par 
Ny 


Re, 

Dans le cas des isotopes pairs, F# = J, et est 
égal à un nombre entier. Les expressions des pro- 
babilités de transition sont de la forme : 

t=2# 2n c\ 4F—2n 
P(F,p,q,t)—= D « (sin 5) (cos à) 


n=1 


en posant 
. æ\2 
(sin 5 = P(1/2, 1/2, — 1/2, 0 
P(1/2, 1/2, — 1/2, t) 


(xH:)? « SL > » 
= — - ; Sin? =V (y)? + (o — w,)? .t 
Date le USE 
Les coefficients af° sont calculés à partir de la 
formule de Majorana. 


Nous pouvons écrire le signal de résonance 


AD7anR foo . œ\?n œ\ 4F—2n ; 
LE se > 1 at (sin?) (cos +) (Na —— N;) Cy er di 


\ 


et après sommation : 


ee) 2n 
s — ES Feu (sin 5) CRTANtE (11) 
T 0 Se 


: < : co ner. 2n 
L'évaluation des intégrales Je ( sin 5) et ds 
0 \ 


est simple, et peut se faire par récurrence ; par 
contre l'évaluation des coefficients X, est fasti- 
dieuse et demande la connaissance des popu- 
lations ÏV, créées par bombardement électronique. 
Ces données étant pratiquement inaccessibles aussi 
bien par le calcul que par l'expérience l’étude pré- 
cise de la forme de raies est impossible. Quelques 
conclusions, cependant, peuvent en être tirées. 
Aux faibles niveaux de radiofréquence, le terme 


JOURNAL DE PHYSIQUE 


en sin? « est le seul important, les autres étant d’un | 


ordre supérieur en yA?, on a donc : 
CPE 
s 4 a. (sin 5 e—t}r dt ; 
1e 0 Z 
À (YH:)° 
1\2 c 
GA) + À) + (e — 0° 


Cette raie est une raie de Lorentz, et sa largeur à ! 


mi-hauteur pour un champ de radiofréquence nul, 
est : 


? 
NO TN = 


Pour une amplitude À, plus grande, les termes 


4 
nc Des : 
en (sin 2) et les termes supérieurs deviennent non 


négligeables, et la raie prend une forme complexe. 
A partir d’une certaine valeur du champ oscillant, 


elle présente un renversement, et l’aspect de la 
courbe est analogue à celui des courbes de réso- 
nance du niveau 6 3#P.. 


En fonction du niveau de radiofréquence, nous ! 


pouvons écrire la demi-largeur sous forme du déve- 
loppement limité : 


AG ASE A2» Pi LE PES 


Dans le cas du niveau 6 3%P,, Brossel a montré 
que les coefficients b ... 
calculé la valeur de a. L’ignorance de la valeur 
des V, ne nous permet pas de préciser ce point. 
Seule l'expérience nous montre, qu’en portant le 


carré de la demi-largeur en fonction du carré du 
champ oscillant, les points sont généralement par- 


faitement alignés. La détermination de Ao, et, par 
suite, de la durée de vie 7, est ainsi très commode. 


Lorsque l’on change la grandeur du champ. 
magnétique oscillant, il se produit une modifi- 
cation de l’intensité lumineuse par suite de la 


variation de la vitesse des électrons et du débit 


électronique. L’étude de l’amplitude à résonance | 


n’est donc pas possible expérimentalement. 


Résultats expérimentaux. — Les trois niveaux - 


identifiés dans le cas du mercure ont été étudiés 


systématiquement, afin de déterminer les facteurs - 
de Landé et les durées de vie. De nombreuses : 


mesures ont été faites à des fréquences différentes 


et nous avons été conduits, pour éviter le risque - 
» = DESERT TE SA . . 
d'erreurs systématiques, à étudier l'influence de la 


structure fine : 
Dans le cas du niveau 6 °F,, le champ de 80 gauss 
nécessaire pour détecter la résonance à 140 Mes 


environ, produit un découplage LS. Les courbes 


représentant les différents états magnétiques du 
niveau 6 %F ne pourront plus être confondues avec 
leurs tangentes à l’origine. On aura donc : 


œ= guB Ho + aH$ + bHS + Me (III) 


les coefficients «4 et b dépendant de la transition 
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sont négligeables, et il a : 
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m — m'. En opérant à fréquence constante en 
faisant varier le champ magnétique,les différentes 
résonances sont décalées et leur superposition pro- 


DéViation du 
galvVanomètre 


[\ 


Déviation du galvanomètre 


> 
72 


Champ : magnétique 
(gauss) 


x A 
68 69 70 


F1G. IL,3. — Résonance d’un niveau #P, du mercure. 
Isotopes pairs, v — 140,4 Mes. 


duira un élargissement apparent de la courbe résul- 


tante. 
—=—_— —— L- r . . 
79,5 80 80,5 g1 Un calcul approximatif [4] montre que, pour que 
& Champ magnétique l'erreur ainsi commise sur la Jargeur limite soit 
(aouss) inférieure à À %, il est nécessaire d’opérer dans un 
F1ic. 11,1. — Résonance du niveau 6%F, du mercure. LEE 
5 Isotopes pairs, v — 140,4 Mcy/s. | Chine 
4 Pour chaque courbe, amplitude du champ oscillant 
7 est donnée en valeur relative. 
È a H°? 
EE 
5 e 
dr. f 28 
s à 
: 2F 17 
À 24 14 
Ps 
ve Fa 10 
44 7 
Pire S 
16 16,50 
Champ magnétique 
(gauss) 
2 
h . 
à = - de Si ; Ha FiG. 11,4. — Résonance du niveau 5d° 6s? 6p ‘D, du mercure 
k Eu CEuR “rs Isotopes pairs, v — 25 Mcys. 


Carré du champ de 
radiofréquence 


Ni 6%F,, v = 25 Mcys. , champ AH, < 40 gauss, soit v < 70 Mc/s. La struc- 

En Ne de a ol largeur à mi-hauteur HE Fes de ee nn 63F te déduite des données 

RC ne D ue numériques de Landolt et Bôrnstein [5]. Aucune 
du champ de radio . 
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influence de cet eflet n’a été constatée pour les 
niveaux %P, et 5d° 6s? 6p 'D,. 

Les figures (11,1 ; 11,3 ; 11,4) montrent quelques 
exemples de réseaux que l’on a obtenus. Sur une 
même figure sont tracées différentes résonances 
correspondant à des niveaux différents de radio- 
fréquence. Le champ oscillant est repéré en valeur 
relative. 

Chaque réseau de,courbe est tracé deux fois 


dans les mêmes conditions, mais en renversant le 


sens du courant magnétisant dans les bobines de | 


Helmholtz. Une moyenne permet ainsi d'éliminer 


le champ magnétique terrestre. 
4 


Une des droites (Aw)? RTE + ayH,)? est 


7) 
représentée sur la figure [1,2.. 


Les différents résultats obtenus peuvent être 


groupés dans le tableau : 


v — 25 Mc/s 
1,248 + 0,003 
(4,49 + 0,15) X 10—7s 


1,121 + 0,001 
(1,58 + 0,16) x 10—68 


NIVEAU v — 140,4 Mc/s 

62F £ 1,2&77 + 0,000% 
d'a 24 (297-042) 10m 

e g 1,4402 + 0,0007 
2 T (2,11 + 0,08) X 1075 

17) £ 1,1203 + 0,0003 
2 T US US ) SCA CES 
Remarques. — 1° La valeur de + trouvée pour le 


niveau 6 °F, dans les expériences à 140 Mes est 
trop faible et est à rejeter : l'erreur provient des 
effets de découplage LS. 

20 D’autres déterminations ont été effectuées 
dans le cas du niveau 6 °F, aux fréquences de 50, 
238, 307,8 et 594 Mc/s. L'accord sur le facteur de 
Landé est toujours excellent à la précision de ces 
mesures, précision qui devient médiocre en champ 
élevé par suite de l’échauffement des bobines de 
Helmholtz produisant une déformation des bobines 
et une modification de l’étalonnage. En outre, la 
mesure de la durée de vie n’a plus aucun sens, l’élar- 
gissement devenant de plus en plus important ; 
à 094 Mc/s la largeur limite est environ 3 fois plus 
grande qu’à 25 Mcys. 

30 Lorsque l’on peut négliger dans l’expres: 
sion (III) les termes en b{H$, l'élargissement par 
effet de découplage est symétrique autour de la 
position déduite de w — gus Hg. La détermi- 
nation du facteur de Landé n’est pas, dans ces 
conditions, affectée par l'élargissement des courbes 
de résonances. 

49 Des mesures ont été faites à des pressions de 
vapeur variées en modifiant la température de la 
goutte de mercure. Aucune variation de la durée 
de vie n’a été constatée (effets de diffusion 
multiple négligeables) [6]. 


IT. Isoropes 1MPAIRS. — L'étude en champ 
faible des différentes résonances des isotopes 199 
et 201 nous a permis, par l’étude des facteurs de 
Landé, de préciser l'identification de certains 
niveaux. En champ plus élevé, le découplage 1J 
n’est pas négligeable, les différentes transitions 
m — m' se produisent à des champs magnétiques 
différents. Chacune de ces raies se traduira par une 
très petite variation du taux de polarisation, le 
rapport signal/bruit est par suite très mauvais, et 


les résonances sont difficiles à identifier au milieu 
du bruit. 


Expérimentalement, l'étude du couplage n’a pu : 


être faite que dans le cas du niveau 6 F4. 

A. Etude du découplage dans le cas de l’iso- 
tope 199 pour le niveau 6 ?F,. 

Résultats expérimentaux. — Ce niveau possède 
deux sous-niveaux hyperfins F = 7/2 et F — 9/2 
dont les facteurs de Landé sont : 


> PSE 1) RUES EM) 


Hs 2F(F +1) | 
1 837 2F(F +1) ; 


En négligeant le deuxième terme, nous pouvons 
écrire : 
RAT 
1,39888. 


Sol2 — 
Eu2 


En champ faible, l’effet de découplage se traduit 
par un élargissement de la résonance, élargissement 
symétrique de part et d'autre du champ magné- 
tique correspondant au facteur de Landé gr. Dans 
ces conditions, nous avons trouvé, à 440 Mc/s : 


gr(9/2) — 1,1082 + 0,0005 gr(7/2) = 1,38850 + 0,0005 
soit 


BE __ 0,8887 + 0,0003. 


niveau 9/2 e 
5 
Valeur théorique : 0,8888 


niveau 7/2 PT — 1,1143 + 0,0008. 
- J M 
Valeur théorique : 1,1111. 


L'accord est donc excellent. 
. À la fréquence de 594 Mecs, les différentes tran- 
sitions m —+ m' sont résolues. La figure II,6 repré- 
sente l’ensemble des résonances du niveau F# — 9/2 
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» 
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dans ces conditions, avec un niveau de radio- 

fréquence élevé. La seule explication compatible 

avec la position des résonances est la suivante : 
Les transitions D sont des résonances doubles 


x 


correspondant à des transitions à deux quanta 
AE 2 

Les transitions 7 sont des résonances Am — 3 
correspondant à la transition à trois quanta. Elles 
n'apparaissent pas sur des courbes tracées à un 
niveau de radiofréquence plus faible. 


_2000 DER Les flèches en pointillé marquent les positions 
ei < déduites pour les résonances simples. Nous n’avons 
-4000 er Eu : pas pu obtenir des courbes où seules les résonances 
- — ue simples apparaissent : lorsque le niveau de radio- 
HE fréquence est faible, le rapport signal/bruit devient 
Pr  e très mauvais, et l’étude devient impossible. 

Te Par une interprétation qualitative des inten- 
de Se sités relatives des trois résonances doubles 
| STE TIREUR 800 1900 9/2 <-+5/2,7/2 <+3/2,5/2 <+ 1/2,et des deuxréso- 
: nances triples, nous sommes conduits à admettre 
Fic. 11,5. — Niveau 6 5F,, Hg. une répartition de population dans les sous-niveaux 

HR 4 4 4 4 OK 4 % 

HO 4 4 VU VU 4 4 2, % 

A : 1 , 

; Ÿ ë ' è ÿ Ÿ ; 

A1 T T T 
nn rt 
NS 9 Mcy 
« : H gau ss 
380 382 384 386 388 390 

_Fic. IL,6. — Résonance du niveau 6 °F,, F = 9/2, Hg, v = 594 Mc/s, D — résonances doubles, T — résonances triples. 


Les flèches en pointillé dans la partie supérieure de la figure indiquent les positions prévues pour les résonances 


simples. 


Les flèches dans la partie inférieure de la figure indiquent les positions théoriques des résonances doubles lors- 


que l’on prend Av = 5 800 Mes. 


DAVIS ES PLATS NS 
nr DE 2 
Fic. 11,7. — En ordonnée, la population relative des 


sous-niveaux Zeeman du niveau 6 2F,, F = 9/2, Hg. 


excités ayant l'aspect de celle donnée par la 
figure [1,7. Une étude précise serait très intéressante 
à ce point de vue, mais les niveaux énergétiques ne 
sont pas, comme dans le cas des isotopes pairs, en 
l’absence du découplage LS équidistants, et écart 
de fréquence entre deux transitions voisines est 
supérieure à leur largeur. Il est donc impossible 
d'employer, pour l’étude de ces phénomènes, les 
formules de Majorana. La méthode de MI Besset, 
Messiah, Horowitz et Winter [7] qui pourrait 
s'appliquer, est dans ce cas d’une extrême com- 
plexité. 

Étude du niveau 9/2: — Les équations des 
courbes composant le diagramme Zeeman sont 
obtenues à l’aide d’un calcul classique [4]. En 
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négligeant le découplage LS, nous sommes conduits 
à diagonaliser l’opérateur. 

D € 

dÙ— - (gg J + £gI 1) Va t- A L.J. 


2mc 


Si l’on pose, en vue de simplifier les équations, 


; 2 él 
Y = v/A et XL — A pee A 
nous obtenons 
He REC) pe VERRE AT 
3 
Me 2 7/2 (y + 2 F4) (y —5 + 3) —2 = 0 (V) 


3 7 
m — + 5/2 (v Het %) (y — 1 F 22) A 0 
Mu ENS 2 (y + 13 
GONE Et 


La construction de ces différentes courbes donne 
le schéma tracé sur la figure I1,5. La structure 
hyperfine Av, écart en champ nul des niveaux 
F —9/2 et F — 7/2, est égale à 4,5 dans-ce sys- 
tème d’unités. Si nous exprimons À en mégacycles, 
nous avONs : 


Av (mégacycles) = 4,5 À. 


Nous avons déterminé expérimentalement l’écart 
entre les résonances — 9/2 <> — 7/2et 9/2 <> 7/2 
du niveau F — 9/2 et à partir des équations précé- 
dentes, par une résolution numérique, nous avons 
évalué Av. L’imprécision des mesures nous conduit 
a 


Av = 5 800 + 100 Mcys. 


Remarque. — Cette détermination déduite des 
équations (V) repose sur l’existence d’une confi- 
guration atomique °F, parfaitement déterminée. 
Le cas du mercure est malheureusement plus com- 
plexe. En outre, nous avons négligé l'influence du 
découplage LS, cet effet s’étant, dans les mêmes 
conditions de champ magnétique, révélé faible dans 
le cas des isotopes pairs. 

L'étude du niveau 7/2 a été également faite de la 
même manière, mais la résolution des différentes 


<— Ft 


200 
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résonances étant moins bonne, l'erreur relative sur 
l'écart des résonances éxtrêmes est plus grande. La 
valeur de Av ainsi trouvée concorde avec celle 
déduite de l’étude du niveau 9/2. 

B. Étude du découplage dans le cas de l'isotope 201 
pour le niveau 6 F4. 

Nous trouvons quatre sous-niveaux hyperfins : 


F = 5/2, 7/2, 9/2, 11/2 


dont les facteurs de Landé théoriques peuvent être 
calculés d’après (IV) : 


11/2 gr = gy X 0,7272. = 0,909 
9/2 gr — gy X 0,8282 — 1,0352 
7/2 gr = gJ X 1,01587 — 1,2697 
5/2 gr — gJ X 1,4285 — 1,7856 


Les structures hyperfines sont plus faibles que 
dans le cas de l’isotope 199, et, par suite, les effets 
de découplage ZJ s’observent dans des champs 
moins élevés. 

Étude expérimentale. — D'une façon générale, 1l 
sera difficile d'étudier les niveaux 7/2 et 9/2. Le 
niveau 7/2 a un facteur de Landé g — 1,267, et sa 


résonance sera superposée à celle des isotopes pairs 


contenus dans l’échantillon étudié. Le niveau 9/2 a 
un facteur de Landé voisin de 1, et les courbes de 
résonance pourront être déformées par des résidus 
de résonance type cyclotron. Des expériences ont 
été faites à 50 Mes, 140 Mes, 238 Mcjs et 
307,8 Mcys : 

— À 50 Mecs, le découplage est très faible, les 
résonances, F—=11/2; PF =9P2 F0 =5)2%$0nt 
observées, et nous obtenons : 


g gr/gJ 


1/2  0,9051 + 0,0005  0,7268 + 0,0005 
9/2 41,0237 + 0,001  0,8280 + 0,0005 
5/2 1,781 + 0,0015  1,4290 + 0,001 


— L'interprétation des spectres obtenus est plus 
complexe que dans le cas de l’isotope 199 (fig. IZ,8). 
Lorsque le découplage est suffisant pour résoudre 
totalement les niveaux 9/2, 7/2 et 5/2, ces différents 
ensembles de résonances se superposent partiel- 
lement, et interprétation devient impossible. Seul 
le niveau 11/2 a pu être étudié à 307,8 Mc/s. Le 
bruit étant important, une étude comparative de 


H Gauss 


F1G. IL,8. — Résonances observées sur un échantillon de mercure enrichi en ?1Hg (88,2 %), niveau 6 3F,, v=— 238 Mes. 
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nombreux enregistrements à dû être faite pour 
déterminer la position des différentes résonances. 
Le diagramme représentant les niveaux éner- 
gétiques en fonction du champ magnétique est 
donné par la figure II,9, et par une méthode ana- 


20 


—+ Fe 
Il 


Q 


4ç0 ROUE EE 00 1000 
2 
Fic. II,9. — Niveau 6 5F,, 201Hp. 


logue à celle employée pour le niveau F — 9/2 de 
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Dans le cas de l’isotope 199, nous avons décelé des 
résonances sur les raies rouges (fig. 11,10) et les 


V =440 Mc4 


Hg 199 
Ps 


Pisotope 199, nous avons déterminé les structures: 


hyperfines : 
Av, (41/2 — 9/2) = 2 860 + 50 Mc/s 
Av; ( 9/2 — 7/2) = 2 340 + 50 Mc/s 
Avs ( 7/2 — 5/2) = 1 820 + 50 Mcys. 


Il est nécessaire de souligner que ces valeurs ont 
été obtenues, — comme dans le cas de l’isotope 199 
—, en négligeant la structure fine. En plus, dans le 
cas de l’isotope 201, nous n’avons pas tenu compte 
du moment quadrupolaire nucléaire. 

Ces différents résultats ne dorvent donc être 
considérés que comme des ordres de grandeurs. 

_ C. Zsotopes impairs. Étude des autres niveaux. — 


F £g EXPÉRIMENTAL 


grlg] EXPÉRIMENTAL 


Hg 199 
isotopes S 
pairs ER 
H gauss 
70 £a 90 100 10 120 
F1G. 11,10. — Résonances observées sur un échantillon de 
mercure enrichi en ‘1%’Hg (73 %), niveau 1D,, 
v = 140 Mes. 
V=ACUAM Ge 
isotopes RËl ie 
Us FE 7 
H gauss 
60 70 80 90 
Fic. 11,11. — Résonances observées sur un échantillon de 


mercure enrichi en 199 Hg (73 %), niveau $ P, v — 140Mc/s. 


raies bleues (fig. I1,11), correspondant respecti- 
vement aux niveaux 1, et 3%P, étudiés précé- 
demment pour les isotopes pairs. Le rapport 
signal/bruit est très faible, et les résultats que nous 
donnons sont par suite assez imprécis. Nousn’avons 
fait, pour la même raison, aucune étude du décou- 
plage IJ. 

Dans le cas des isotopes 201, nous n’avons décelé 
aucune résonance. Les différents résultats sont 
réunis dans le tableau : | 


gr[gJ THÉORIQUE 


3/2 
"Pa 5/2 1,15 + 0,0005 0,7995 4 0,0008 0,8 
ip, 32 1,3 + 0,0 1198 + 0,01 1,2 
2 5/2 0.902 + 0,005 0805 æ+ 0.01 08 


Expériences sur le sodium. — Le queusot d’un 
tube du type n° 2 est rempli de sodium soigneu- 
sement distillé. Ce tube, ainsi que les bobines pro- 
duisant le champ de radiofréquence, sont placés 
dans un four thermostaté. La température opti- 
mum est d'environ 1800. À des températures plus 
élevées, le verre de l'enceinte est rapidement 
attaqué. 

Nous avons étudié différentes raies du spectre 
du sodium non absorbées par le verre du tube et 
dont les longueurs d’onde correspondent au 
domaine de sensibilité des photomultiplicateurs. 


Un taux de polarisation appréciable n’a été observé 
que sur les raies 52D;y, — 32Py2 (4983 À), 
5 2Djo — 3 2P,y2 (4 979 À) et 4 Dj — 8 2Pajo 
(5 688 À), et nous avons décelé les résonances des 
niveaux 5 2D,2, 4 2Dip et 5 2D.,, (fig. 11,12). 

Des expériences préliminaires ont été faites 
à 25 Mc/s. Les facteurs de Landé, déterminés à 
partir de la position de ces résonances, nous 
montrent que le découplage IJ est à peu de chose 
près total. Les structures hyperfines de ces niveaux 
sont donc faibles par rapport à la fréquence de 
25 Mcys. 
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Différentes expériences ont été faites à 50 Mcys. 
Nous constatons que la largeur limite en champ de 
radiofréquence nul est plus grande qu'à 25 Mecys. 
Cet effet est dû au découplage ZS, l'écart des 
niveaux Do — Days étant très faible (distance des 


niveaux D 2), — 52),,41:= 0.02%/cmvdis; 
tance 4 2D,,, —42D,, —=0,083cm°*. 
#Ds,, 
22 20 18 16 14 12 Gauss 


5? D, 
on 
me Se 


1 —_—) le 1 a — 
21; 20 16 12 Lauss 
Fc. 11,12. — Résonances observées sur le sodium. 
V2 NIC/S: 


Courbe 1, en haut : raies 5688 et 5683 À. 
Courbe 2, en bas : raies 4983 et 4979 À. 


La détermination de la durée de vie est donc 
complexe et nécessite une étude complète de leffet 
Paschen-Back hyperfin et fin. Ce travail est actuel- 
lement en cours et les résultats seront donnés dans 
une publication ultérieure. 


La méthode d’excitation électronique et ses 
possibilités. — Sur le mercure et sur le sodium 
nous n'avons décelé et étudié la résonance magné- 
tique que de quelques niveaux seulement. Nous 
pouvons en préciser les raisons et définir ainsi le 
domaine d’application de la méthode d’excitation 
électronique : 


1) AMPLITUDE DU CHAMP DE RADIOFRÉQUENCE. 
—- Pour que des transitions entre deux sous-niveaux 
magnétiques puissent être observées, il est néces- 
saire que la valeur de la probabilité de transition 
soit appréciable au bout d’un temps de l’ordre 
de la durée de vie +. La valeur numérique du 
champ À, nécessaire peut être calculée à partir des 
formules de Majorana. On trouve ainsi une ampli- 
tude de l’ordre du gauss pour des durées de vie 
de 10-7 sec et d’une dizaine de gauss pour des 
durées de vie de 10-$ sec. Dans ce dernier cas, les 
trajectoires électroniques sont violemment per- 
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turbées, le taux de polarisation des raies émises est 
par suite très faible et il devient difficile d'empêcher 
les décharges dans l’enceinte. 

La méthode d’excitation électronique, dans le 
cas des transitions dipolaires magnétiques, ne peut 
par suite s'appliquer qu'aux niveaux dont la duree de 
pie est superieure à, approzimativement, 5.10 
seconde. 

Dans l’atome de mercure, la plupart des niveaux : 
ont une durée de vie très courte. Pour le 
niveau 61P,, on a r <10-# seconde, pour le 
niveau 7 %$,, r # 108 seconde [8], et il en est de 
même des niveaux D [9]. Il est aisé de comprendre 
pourquoi peu de résonances ont été observées. 


2) CRÉATION D’INÉGALITÉ DE POPULATION DANS 
LES ÉTATS EXCITÉS. — Dans la première partie, 
nous avons signalé qu’il était impossible de donner 
des règles générales permettant de prévoir, même 
de façon qualitative, les taux de polarisation. Il 
semble toutefois que, lorsqu'un choc électronique 
porte un atome de l’état fondamental, de multi- 
plicité F, à un état excité, de multiplicité F” plus 
élevée (F' > F +2), il y a création d’inégalité de | 
population dans cet état excité, le choc électro- 
nique défavorisant une excitation se faisant avec 
APR EE MO 

Les niveaux que nous avons étudiés (%%,, %P», 
1D, de l'atome de mercure, et D;,, de l'atome de 
sodium) satisfont cette condition. 


Conclusions. — Ce travail d'exploration a donc 
fourni des résultats positifs et il nous montre que la 
méthode d’excitation électronique complète, pour 
l'étude des niveaux élevés, la méthode d’excitation 
optique, de beaucoup préférable lorsqu'elle peut 
s'appliquer. 

Nous pouvons envisager diverses extensions : 

1) Dans le cas d’un atome à spin nucléaire, il 
doit être possible de détecter les résonance AF —1. 
Un montage permettant d'étudier ces transitions 
pour le niveau 6 %F, de lisotope 201 du mercure 
est en cours de réalisation. Nous pourrons ainsi 
déterminer de façon précise les structures hyper- 
fines dont nous avons donné un ordre de grandeur. 
Les résultats de cette étude”seront donnés ulté- 
rieurement. 

2) Par une technique tout à fait analogue à 
celle décrite ci-dessus, nous pourrons étudier, en 
bombardant une vapeur, un certain nombre 
d’atomes : potassium, césium, rubidium. Le cas des 
gaz rares pourra être également envisagé. Le pro- 
blème, cependant, sera plus complexe, les spectres 
comportant un très grand nombre de raies. 

D’autres éléments plus difficilement volatils 
pourront être étudiés par une technique de jet 
atomique. Afin d’avoir un flux lumineux compa- 
rable à celui dont on disposait avec une vapeur, il 
sera nécessaire d'opérer sur des jets atomiques de 
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grande densité. Il faudra mettre au point un mon- 


_tage à pompage continu en prenant des précautions 


particulières pour ne pas contaminer la cathode (*). 

3) Les expériences que nous avons faites sur le 
sodium montrent que, pour cet élément, peu de 
raies émises par bombardement électronique sont 
polarisées. Si l’on crée une inégalité de population 
par pompage optique dans le niveau fondamental, 
nous pouvons penser qu’elle sera conservée lorsque 
les atomes seront portés dans un état excité par 
bombardement électronique. On peut espérer ainsi 
étudier de nombreux niveaux. 

4) Par pompage optique, il est possible d’enri- 
chir sélectivement les sous-niveaux magnétiques de 
l’état fondamental. Un effet analogue peut être 
obtenu par bombardement électronique. En effet, 
si nous réglons la tension d'accélération des élec- 
trons de façon à n’exciter dans des conditions de 

(*) Un tel montage a été utilisé par H. Bruck [10]. 
Notons toutefois que la cathode qu’il employait, était 
moins fragile que les cathode industrielles, au rendement 
beaucoup plus élevé, que nous utilisons. 

E. Geneux et Mlle B. Vincenz [11] ont, par cette tech- 


nique, récemment mis en évidence la résonance magné- 
tique de plusieurs niveaux excités du cadmium. 
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polarisation favorable que le premier niveau excité, 
il est vraisemblable que, lorsque les atomes excités 
retomberont à l’état fondamental, une inégalité de 
population sera créée dans cet état. 

Le bombardement électronique peuple de façon 
symétrique les niveaux —m et + m, aussi, 
contrairement à l’excitation optique en lumière 
polarisée circulaire qui crée une orientation, nous 
ne pourrons obtenir qu’un alignement pour les 
niveaux fondamentaux tel que F > 1. 

Pour obtenir un alignement appréciable, il est 
nécessaire, d’une part, que les raies qui corres- 
pondent à une transition entre un niveau excité et 
le niveau fondamental soient polarisées — le 
sodium est donc exclu — et, d’autre part, que les 
collisions sur les parois et sur les électrodes ne 
détruisent pas l’alignement, l’excitation électro- 
nique ne permettant pas l'emploi d’un gaz tampon. 
Cette dernière condition n’est pas remplie pour les 
métaux alcalins ; le cas de l’isotope 201 du mercure 
semble se présenter plus favorablement [12], si les 
chocs sur les parois métalliques des électrodes ne 
détruisent pas l’alignement. 

Manuscrit reçu le 25 mars 1959. 
Extrait de [13]. 
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SUR LES SYSTÈMES DE CORPUSCULES 
EN THÉORIE FONCTIONNELLE NON RELATIVISTE 


Par FLorence AESCHLIMANN, 


Théories physiques, Institut Henri-Poincaré, Paris. 


Résumé. -- Examen des conditions générales conduisant à une théorie fonctionnelle non rela- 
tiviste des systèmes de corpuscules. Représentation de chaque corpuscule d’un système par une 
fonction u;(P, T) de l’espace et du temps. Ordination formelle des corpuscules d’un système tenant 
compte à la fois de l’indiscernabilité des corpuscules de même espèce et de la nécessité d'exprimer 
le principe d’exclusion pour les fermions. Utilisation de la théorie générale des prévisions ; déter- 
mination des arguments de la fonctionnelle des prévisions, introduction de l’espace fonctionnel de 
configuration, détermination de la forme générale de l'équation prévisionnelle, étude des conditions 
de raccordement avec la mécanique ondulatoire usuelle des systèmes et avec le cas d’un corpuscule 
unique en théorie fonctionnelle. Énoncé du principe fondamental des équations des corpuscules ; 
chaque fonction w,(P, T) représentant un corpuscule du système considéré obéit à une équation 
non linéaire dépendant des autres fonctions w,{P, T) du système. Étude sommaire des conditions 
de raccordement avec la mécanique classique. 


Abstract. — À study of the general conditions required to set up a non-relativistical functional 
theory for systems of particles. Representation of each particle of a system by a space-time func- 
tion u;(P, T). Formal classification for the particles of the system with respect to indiscernability 
and the exclusion principle for the fermions. Using the general theory of predictions, the argu- 
ments of the functional for predictions are fixed, functional configuration-space is introduced and 
the general form of the equation for the predictional elements is determined. Conditions incurred 
by the limiting cases of the usual wave mechanics for systems and of the single particle in the 
functional theory of particles. Fundamental principle for the equations of motion concerning 
the particles of a system : each function w;(P, T) describing a particle of the system obeys a 
non-linear equation which depends upon the functions w(P, T) of the other particles of the 
system. Brief investigation of the limiting case of classical mechanics. 


1. Introduction. — L'étude des systèmes de cor- 
puscules en théorie fonctionnelle non linéaire pose 
des problèmes très difficiles qui n’ont pas encore été 
abordés. Le cas de la théorie non relativiste est 
évidemment celui par lequel il convient d’essayer 
de commencer cette étude car, d’une part il permet 
de laisser de côté toutes les questions de propa- 
gation des interactions et, d’autre part, il permet 
de considérer un système d’un nombre déterminé 
de corpuscules. Il est bien clair que l’on ne peut se 
limiter à l’étude d’un corpuscule unique dans un 
champ donné et que tout problème ayant un inté- 
rêt physique fait intervenir un système de corpus- 
cules. En théorie des systèmes apparaissent des 
circonstances qui ne se présentent pas dans le cas 
d’un seul corpuscule et c’est là ce qui fait l’intérêt 
de l’étude des systèmes, même en théorie non rela- 
tiviste. 


2. La théorie fonctionnelle des corpuseules. — 
En théorie fonctionnelle [1], un corpuseule C est 
représenté par une fonction de point de l’espace 
physique et du temps, soit u(P, T) dite onde phy- 
sique, tandis qu’en mécanique classique comme en 
mécanique ondulatoire usuelle, un corpuscule est 
représenté par un point M de l’espace physique 
qui est supposé euclidien à trois dimensions. Cette 


hypothèse fondamentale de la théorie fonctionnelle 
permet de traiter le cas d’un corpuscule dans un 
champ donné ; mais comment généraliser cette 
hypothèse dans le cas d’un système ? C’est Jà une 
première question fondamentale préalable. 

On admet de plus que la fonction u(P, t) appar- 
tient à un espace fonctionnel séparable (R,) [2] et 
qu’elle obéit à une équation non linéaire 


fu = Q(u, v) (2) 


où Q est un terme non linéaire en w et ses dérivées 
et où est un opérateur linéaire tel que 


F#V—=0 (2) 


est l’équation d’onde de la mécanique ondulatoire 
usuelle du problème correspondant. Si alors on 
considère un système, comment ces conditions 
seront-ellesg énéralisées ? Nous avons laune seconde 
question fondamentale préalable. Examinons 
d’abord ces deux questions. 


3. Comparaison avec la théorie de Schrôdinger. 
— En mécanique ondulatoire usuelle dans le cas 
d’un système, d’après la théorie développée par 
Schrüdinger [3], il y a pour le système une seule 
onde Ÿ qui se propage dans l’espace de configu- 


# 
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ration du système en obéissant à une équation 
d'onde linéaire 


Pie 0: 


La condition que nous avons énoncée ci-dessus 
concernant l’équation de l’onde physique w dans 
le cas d’un corpuscule unique ne peut subsister 
inchangée dans le cas d’un système : ou bien la 
liaison entre l’équation (1) en et l’équation (2) 
en % exprimée par le même opérateur £ disparait, 
ou bien on aura pour le système une seule onde x 
qui se propagera dans l’espace de configuration. 

Comme l’a fait remarquer depuis longtemps 
M. Louis de Broglie [4], l’espace de configuration 
est un espace fictif, une construction mathématique 
auxiliaire que la mécanique classique a introduit 
pour l’étude d’un système de points matériels en 
mouvement dans l’espace physique. De même en 
microphysique, c’est dans l’espace physique qu’il 
faut considérer les corpuscules, quelle que soit la 
façon dont on les représente. 

En mécanique ondulatoire usuelle, un corpus- 
cule est figuré par un point variable hbre dans 
l’espace physique (point M dont les coordonnées 
sont 3 variables libres x, y, z) ; pour un système 
de n corpuscules dont les n points M; figurent les 
corpusCules, on peut former un point variable 
libre M de l’espace de configuration. La fonction 
d’onde % a une signification purement prévi- 
sionnelle ; à partir de cette fonction, on calcule 
toutes les prévisions concernant le système. On doit 
donc avoir une seule fonction d pour tout le sys- 
tème (et non une fonction d; pour chaque corpus- 
cule C;) car une évolution dans le temps n’a de 
sens que pour la totalité du système et non pour 
chaque partie le constituant. Les arguments de 
cette fonction 4 doivent être les points figurant les 
corpuscules du système, soit 


b—= M, M, ..., Mat). 


On peut considérer la fonction 4 comme une 
onde se propageant dans l’espace de configuration. 
D'autre part, en mécanique classique, le mouve- 
ment d’un système est régi par une seule fonction, 
la fonction de Jacobi S de l’espace de configuration; 
comme à l’approximation de l’optique géomé- 
trique, la description de la mécanique ondulatoire 
doit se raccorder avec la description classique, il 
faut que le système soit décrit par une seule fonc- 
tion Ÿ se propageant dans l’espace de configuration. 
Ceci vient justifier le fait de considérer l'onde Ÿ 
dans l’espace de configuration. 


4. Le premier principe de la théorie fonction- 
nelle, — En théorie fonctionnelle, on admet tout 
d’abord qu’un corpuscule C est représenté par une 
onde physique u(P, T) (de même qu’en mécanique 
ondulatoire usuelle, il est représenté par un 
point M), qui se propage dans l’espace physique. 
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Un corpuscule n’est jamais complètement séparé 
des autres corpuscules ; au contraire la notion de 
système a quelque chose d’artificiel, elle n’apparait 
que par une coupure établie entre les corpuseules 
du système et le reste de l’univers. Donc on ne 
peut changer qualitativement la représentation 
d’un corpuscule, qu’on le considère comme isolé, 
comme appartenant à un système déterminé, ou 
comme appartenant à tout l’univers : si dans l’un 
des cas, un corpuscule est représenté par un. 
point M, il demeurera représenté par ce point M 
dans tous les cas ; si c’est une fonction u(P, T) 
qui le représente dans l’un des cas, c’est une telle 
fonction qui le représentera dans tous les cas. Nous 
sommes ainsi conduits à admettre ce premier prin- 
cipe fondamental : 


PRINCIPE DE LA REPRÉSENTATION : Les corpus- 
cules C d’un système physique S seront représentés 
chacun par une fonction u(P, T). 

Ainsi se trouve donnée une réponse à la première 
question préalable. 


9. Ordination des corpuseules. — La mécanique 
ondulatoire usuelle des systèmes nous a appris que 
les corpuscules de même espèce sont indiscernables 
et qu’une ordination de corpuscules de même espèce 
est dénuée de signification physique. Cependant, en 
mécanique ondulatoire on numérote les corpuscules 
et les points variables libres qui les représentent, 
sachant que les quantités ayant une signification 
physique seront définies à une permutation des 
corpuscules de même espèce près, donc à une per- 
mutation près des points variables qui les repré- 
sentent. En effet, les points variables représentant 
les corpuscules figurent comme arguments dans les 
fonctions d'ondes 4 du système. Un numérotage se 
trouve fixé par les places des arguments dans la 
fonction : il y a nr places d'arguments pour des 
points variables dans une fonction à du système. 
Le point variable M mis à la je place prend de ce 
fait un numéro par la place qu’il occupe, on le 
notera M; ; il doit être distinct du point variable 
qui occupe une autre place, par exemple la Æ® 
place ; même s’il représente un corpuscule de même 
espèce que M,, il est un point variable différent 
qu’on notera M,. Cette ordination formelle est indis- 
pensable à la représentation mathématique du 
système, mais cela ne lui confère aucune signi- 
fication physique ; il y a un ordre dans les points 
variables figurant les corpuscules comme:il y a un 
ordre dans les mots d’un texte, cet ordre ne se 
transmet pas dans les choses mais il est indispen- 
sable à la description théorique. On respectera 
l’indiscernabilité des corpuscules de même espèce 
en fixant que les propriétés physiques seront défi- 
nies à une permutation près des corpuscules. 

De plus l’ordination des points figurant les cor- 
puscules est indispensable pour la formulation de 
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certaines lois physiques, par exemple pour expri- 
mer le principe d'exclusion de Pauli lorsqu'on pose 
que la fonction d’onde 4 change de signe-si Pon 
permute deux fermions. Pour que l’on puisse per- 
muter, il faut que la permutation soit formulable, 
qu’on ait préalablement fixé formellement un ordre 
permettant de distinguer le cas où la permutation 
est effectuée de celui où elle n’a pas été effectuée, 
de façon qu’on sache quel signe mettre devant la 
fonction d. S'il n’y avait pas ordination formelle 
des corpuscules, une permutation des arguments 
de Ÿ ne serait pas descriptible, il ne pourrait pas y 
avoir changement de signe de d lors d’une permu- 
tation, la fonction d’onde 4 serait symétrique. 
Ainsi une ordination formelle des éléments repré- 
sentant les corpuscules est essentielle pour la formu- 
lation du principe d'exclusion de Pauli. 


6. Prévisions en théorie fonctionnelle. — Dans 
les raisonnements précédents, le fait que les cor- 
puscules soient représentés par des points variables 
libres en mécanique ondulatoire usuelle n’inter- 
vient pas d’une façon effective. Les raisonnements 
doivent être répétés sans changement si l’on repré- 
sente les corpuscules par des fonctions u(P, T), 
mais ce qui intervient essentiellement, c’est que 
Jon peut calculer des prévisions sur le système. 
C’est sur les fonctions d’ondes d servant à calculer 
les prévisions que se formulent les conditions de 
symétrie ou d’antisymétrie (formulation du prin- 
cipe d’exclusion) en mécanique ondulatoire, non 
sur des éléments intrinsèques au système. 

Le but minimum que l’on peut fixer à une théorie 
physique est de permettre le calcul de prévisions 
(même si l’on exige davantage d’une théorie phy- 
sique) car s’il n’y a pas prévisions, il ne peut y 
avoir vérification expérimentale et par suite con- 
nexion avec l’expérience physique. Mais ce but 
minimum est quand même le but pratique le plus 
important, même si l’on se fixe comme but fonda- 
mental une explication rationnelle, si bien que dans 
toute théorie physique il y aura des éléments per- 
mettant le calcul de prévisions ; la théorie fonc- 
tionnelle des corpuscules ne peut échapper à cette 
condition. D’autre part, de même que la méca- 
nique ondulatoire est tenue de se raccorder à la 
mécanique classique au moyen de l’approximation 
d'optique géométrique, la théorie fonctionnelle doit 
se raccorder à la mécanique ondulatoire usuelle par 
l’approximation qui consiste à remplacer la repré- 
sentation fonctionnelle par la représentation ponc- 
tuelle. Pour formuler ce raccordement de façon 
précise, 1l suffit d'utiliser les. résultats de la théorie 
générale des prévisions : à chaque observation 
maximale du système S considéré, effectuée à un 
instant t,, on peut faire correspondre un élément 
initial X,, et à cet X, correspond à l'instant t{ un 
élément de prévision X(1). À l’approximation de la 
mécanique ondulatoire, X, se réduit à une” fonc- 


No 8-9 
tion d’onde initiale ‘4, et X(t) à une fonction 
d'onde v(t) du système. Un élément de pré- 
vision X(t) dépend du système S considéré, c’est-à- 
dire dépend de l’ensemble des 7 corpuscules cons- 
tituant le système S. Un élément X(1) est un être 
mathématique ; s’il dépend des » corpuscules du 
système, ce ne peut être que par l'intermédiaire de 
leur représentation mathématique ; c’est donc une 
fonction de n éléments représentant les n corpus- 
cules du système S. Comme nous avons adopté une 
représentation fonctionnelle, X(t) doi être une 
fonctionnelle des n fonctions u,(P, T), u{(P, T), 

., a(P, T) représentant les n corpuscules du sys- 
tème S, Soit 


xt = RAP, The) D), halle Ti 


Les conditions de symétrie (bosons) et d’anti- 
symétrie (fermions) s’exprimeront [5] sur les élé- 
ments X exactement comme en mécanique ondu- 
latoire usuelle sur les fonctions 


dCW;, M, PAC) M,, t). 


7. Espace fonctionnel de configuration. — L’es- 
pace de configuration de la mécanique ondulatoire 
usuelle est formé par le produit cartésien 
de 7 espaces euclidiens à 3 dimensions, chacun. 
correspondant à l’un des corpuscules, c’est-à-dire 
à l’un des points variables représentant un cor- 
puscule du système, chacun de ces points variables 
parcourant l’espace euclidien pris comme espace 
physique. 

En adoptant la notation de la somme directe ©, 
le point variable M représentant le système dans 
l’espace de configuration sera donc formé à partir 
des points M; représentant les corpuscules par la 
relation 


M=M, OM, © .… OM; ® .…. ® M. 


Ce qui est constitué par un point variable de 
l’espace de configuration, c’est donc l’ensemble des 
éléments mathématiques représentant les cor- 
puscules du système en mécanique ondulatoire. 

En théorie fonctionnelle, les éléments mathé- 
matiques représentant les corpuscules du système 
sont les r fonctions u;(P, T), par suite”c’est l’en- 
semble de ces éléments qui doit constituer l’espace 
de configuration du système. Chacune des fonc- 
tions u;(P, T) appartient à l’espace (R,) des fonc- 
tions u. l’espace de configuration (C) sera donc ici 
un espace fonctionnel, et le point u de cet espace 
représentant le système S sera la somme des n fonc- 
tions u(P, T) représentant les corpuscules, soit : 


; CISAR TERRES (Ra) 
u(P,T) = u(P,T) QŒu,(P,T) 
(>) ..® u;(P, T) (@) 5,000 ® ua(P,T). 


Alors, de même qu’en mécanique ondulatoire 
usuelle, la fonction d'onde Ÿ du système a pour 
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arguments le point M et le temps t de l'horloge de 
l’observateur, soit 4(M, t), en théorie fonctionnelle, 
l'élément de prévision X aura pour arguments le 
point w de l’espace fonctionnel de configuration (C) 
et le temps t de l'horloge de l'observateur, soit 
X[u(P, T), t]. 


8. Les équations prévisionnelles. — I1 nous faut 
maintenant chercher une réponse à la seconde 
question préalable et déterminer des équations 
pour l’évolution du système. Nous ne sommes 
guère guidés dans cette recherche que par des 
conditions de raccordement : d’abord, raccor- 
dement avec la mécanique ondulatoire usuelle des 
systèmes, ensuite avec la théorie fonctionnelle du 
corpuscule unique, enfin avec la mécanique clas- 
sique des systèmes. 

Dans une théorie nouvelle, il reste toujours des 
traces de la structure de la théorie plus ancienne ; 
c’est en cela que consiste la persistance des struc- 
tures formelles [6]. En théorie fonctionnelle des 
systèmes de corpuscules, il restera donc des élé- 
ments de la structure formelle de la mécanique 
ondulatoire des systèmes usuelle. Mais cette méca- 
nique ondulatoire est purement prévisionnelle ; ce 
n’est donc que sur le plan prévisionnel que cette 
mécanique peut nous fournir des conditions et de 
telles conditions seront très insuffisantes pour par- 
venir à une réponse à la question posée. Examinons- 
les cependant tout d’abord. 

Les éléments de prévision X({) sont la généra- 
lisation pour la théorie fonctionnelle des fonctions 
d'onde d de la mécanique ondulatoire usuelle ; 
c’est donc sur ces éléments X(t) que se portera ce 
qui subsiste en théorie fonctionnelle de la structure 
de la mécanique ondulatoire. Tout d’abord, nous 
pouvons faire intervenir les conditions qui résul- 
tent de la théorie générale des prévisions. En pre- 
mier lieu, les éléments de prévision X({) appar- 
tiennent à un espace abstrait (V) qui en général 
n’est pas un espace de Hilbert, mais qui contient 
au moins un espace de Hilbert. D’autre part, nous 
pouvons introduire un opérateur d'évolution 
‘U(4, to) tel que 


X(t) = Us, to) Xo ; (3) 


. toute équation d'évolution des prévisions peut alors 


s’exprimer par une équation sur l'opérateur d’évolu- 
tion U(é, t,). En mécanique ondulatoire des sys- 
tèmes usuelle, l'équation d'évolution dérive immé- 
diatement de l’équation d’ondes ; on a 


., DU 
Ho Umo = -— D (4) 


où A est l’opérateur hamiltonien du système, soit 


An pu Mat). 
ji 2j 


SUR LES SYSTÈMES DE CORPUSCULES 733 


Où Vuo( Mi, Mo, ..., M,, t) est le potentiel de la 
mécanique classique et où 


Pi =a ti V; 


s’il n’y a pas de champ magnétique, et où 
. £ 
Pi ai V; —"A 


en présence d’un champ magnétique dérivant d’un 
potentiel vecteur A. En désignant par Two l’opé- 
rateur énergie cinétique du 7° corpuscule on peut 
encore écrire 


Huo = 2 Tuos + V. 
3 


La théorie générale des prévisions conduit 
pour ‘U à l’équation intégrale [7] : 


t 
As, to) = UE, to) + JF, 5) Ur, to) de 


où ‘U0(4, t,) est l'opérateur d'évolution d’un sys- 
tème fictif, dit substratum de S, qui obéit à une 
équation de la forme 
; ., DU 
Ho ON Sa 
En théorie non relativiste des systèmes, sans 
création ni annihilation de corpuscules dans le 
système S étudié, on doit admettre que l’opéra- 
teur ‘U est dérivable et que par suite # est déri- 
vable, et ceci conduit à l'équation différentielle 
AU 
(H, + R(t)) WU = — ne (5) 


avec F 
RIDE) 
En comparant avec l'équation (4) de la méca- 
nique ondulatoire, on voit que l’on est conduit à 
poser 


H = H, + R(t) 


et à admettre que l’opérateur U en théorie fonc- 
tionnelle obéit à une équation de la forme (5) 
Comme opérateur /1,, on peut prendre celui corres- 
pondant au substratum S, formé par les mêmes 
corpuscules que S, mais qui sont supposés sans 
interaction entre eux ; ceci conduit à poser 


Ho = À Ti. 
j 


Au lieu de considérer l’équation (5) entre opéra- 
teurs, en utilisant la relation (3), on obtient une 
équation sur les éléments de prévision À. Ainsi 
finalement, la persistance des structures de la 
mécanique ondulatoire usuelle en théorie fonc- 
tionnelle nous conduit au principe suivant : 


PRINCIPE DES PRÉVISIONS : Pour un système S de 
corpuscules en théorie fonctionnelle non relativiste 
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les fonctionnelles de prévision X(u, t) obéissent à une 
équation de la forme 


dX 


HX = — 1à 


avec 
H=Z2T;+R 
e) 


où T; serait l'opérateur H correspondant au cas 
du j® corpuscule considéré isolément et où R est un 
terme lié à l'interaction des corpuscules entre eut. 

La fonctionnelle X a pour arguments un élément u 
de l’espace fonctionnel de configuration (C) et le 
temps t de l'horloge de l'observateur. Cette fonc- 
tionnelle X est un élément de l’espace (VU) des élé- 
ments de prévision. L'espace (VU) contient au moins 
un espace de Hilbert. Les opérateurs T;, R, H sont 
des opérateurs de l’espace (V). 


9. Les équations fondamentales. — Le principe 
précédent exprime ce qui subsiste en théorie fonc- 
tionnelle de la structure de la mécanique ondula- 
toire usuelle des systèmes, mais ces conditions pré- 
visionnelles ne contiennent aucun élément propre 
à la théorie fonctionnelle des corpuscules ; le fait 
qu’un corpuscule soit représenté par une fonction 
n'intervient pas d’une manière essentielle. Dans le 
cas d’un corpuscule unique (soumis ou non à un 
champ extérieur), nous avons l’équation (1) à 
laquelle doit obéir toute fonction u représentant le 
corpuscule. Cette fonction w dépend du corpuscule 
et des circonstances dans lesquelles il est placé, en 
d’autres termes la fonction # dépend du corpuscule 
et du reste de l’univers, ce reste de l’univers mani- 
festant son influence par le terme V contenu dans 
l'opérateur #, par le potentiel vecteur À et par le 
terme non linéaire Q. Si au lieu de considérer le 
corpuscule comme étant isolé, on le considère 
comme un élément d’un système S, ceci ne change 
aucune condition physique, mais change seulement 
le sens de ce que nous appelons système et de ce que 
nous appelons le reste de l'univers. Il doit donc 
obéir encore à la même équation, mais dans les 
termes V, À, Q, nous devons distinguer une partie 
due à l’influence de l’extérieur du système S et une 
partie due à l’influence des autres corpuscules du 
système. Ceci naturellement concerne le cas d’inter- 
actions assez faibles pour que l’on puisse, avec une 
adéquation suffisante, considérer un corpuscule 
comme isolé. Le même raisonnement de raccor- 
dement avec le cas d’un corpuscule isolé peut être 
appliqué à chaque corpuscule du système. Mais si 
les interactions se renforcent, il n’y a aucun élé- 
ment qui puisse venir détruire l’équation 


Fj uj = Q; 


du 7° corpuscule du système ; seule la forme des 
termes peut changer. Aïnsi les conditions de raccor- 
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dement avec le cas d’un corpuscule unique conduit 
à admettre ce principe fondamental : 


PRINCIPE DES ÉQUATIONS DES CORPUSCULES : 
Chaque corpuscule C; d'un système S est représenté 
par une fonction u(P, T) qui obéit à une équation 

Fu = Q,; (6) 


où £, est un opérateur linéaire de l’espace (R,) et Q; 
un terme non linéaire qui est un élément de l’es- 
pace (R,). L'opérateur #; et le terme Q; dépendent 
des fonctions u,(P, T) représentant les autres cor- 
puscules du système. 

Ainsi on a un système de 7 équations en u,(P, T), 
u{(P, T), ..., u(P, T) pour caractériser l’évolu- 
tion d’un système de 7 corpuscules. Aux raisons 
générales données ci-dessus concernant l’adoption 
de ce principe viennent s'ajouter les raisons données 
par M. Louis de Broglie dans le cadre plus strict 
de la théorie de la double solution [8]. 


10. Raccordement avec la mécanique classique. 
— La théorie fonctionnelle des systèmes de cor- 
puseules doit aussi se raccorder avec la mécanique + 
classique. Le raccordement concernant l’aspect pré- 
visionnel est automatiquement satisfait par le 
raccordement avec la mécanique ondulatoire 
usuelle, puisque cette mécanique se raccorde par 
approximation d'optique géométrique avec la 
mécanique classique. Mais par le principe précé- 
dent, nous venons d'admettre un postulat beau- 
coup plus fort que ceux de la mécanique ondu- 
latoire usuelle. Il faut donc examiner s’il y a bien 
raccordement avec la mécanique classique. Le prin- 
cipe précédent fixe qu’il y a une équation pour 
chaque corpuscule du système. Le raccordement 
s'effectuera d’abord en remplaçant la représen- 
tation fonctionnelle par une représentation ponc- 
tuelle. Dans ce cas, les équations (6) vont se trouver 
remplacées par des équations sur les points M, 
représentant les corpuscules. Nous obtiendrons 
donc un système d'équations correspondant à 
chaque corpuscule ; ceci est bien en accord avec le 
principe des équations différentielles de la méca- 
nique classique des systèmes de points matériels ; 
ces équations s’explicitent sous la forme 

> — 
m L’; = F;. (7) 

La condition de raccordement avec la mécanique 
classique exige donc qu’à partir des équations (6), 
en effectuant les approximations convenables, on. 
retombe sur les équations précédentes. Ceci vient 
fournir des conditions sur les opérateurs £ .; nous 
les examinerons ultérieurement. Mais on voit déjà 
que la condition de raccordement avec la méca- 
nique classique fournit un argument de plus pour 
qu'il y ait une équation pour chaque corpuscule, 
d’où un système d’équations (6) pour un système 
de corpuscules, et non pas une seule équation. 
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En mécanique ondulatoire usuelle, on doit aussi 
avoir une persistance de la structure introduite par 
les équations (7) ; mais comme cette théorie est 
purement prévisionnelle, avec représentation ponc- 
tuelle des corpuscules, les équations (7) doivent 
subsister comme équations entre opérateurs ; c’est 
bien ce que l’on constate [8]. Ici, par la repré- 
sentation fonctionnelle des corpuscules, il apparaît 
un aspect de la théorie qui n’est pas prévisionnel 
et qui exige un raccordement de nature différente, 
c’est pourquoi le raccordement avec la mécanique 
ondulatoire usuelle ne suffit pas et que la méca- 
nique classique vient imposer directement d’autres 
conditions. 

Il ya enfin d’autres conditions de raccordement, 
ce sont celles avec Ia théorie de la double solution à 
l’approximation linéaire, théorie étudiée par MM. 
Louis de Broglie et Andrade e Silva [9] ; nous les 
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étudierons après avoir terminé l'étude des condi- 
tions de raccordement avec la mécanique classique. 


Conclusion. — La conclusion de cette étude est 
donc qu’on peut construire une théorie fonction- 
nelle non relatiste des systèmes de corpuscules qui 
satisfasse aux conditions de raccordement avec la 
mécanique ondulatoire usuelle des systèmes et avec 
la mécanique classique. Dans cette théorie, chaque 
corpuscule C; du système étudié est représenté par 
une fonction u;(P, T) de l’espace et du temps appar- 
tenant à un espace fonctionnel séparable (2,). 
Chaque fonction u;(P, T) obéit à une équation de la 
forme (6). Un système de » corpuscules a donc son 
évolution régie par un système de » équations de la 
forme (6). La forme précise de ces équations est 
donnée à partir des conditions de raccordement. 
avec la mécanique classique et avec le cas du cor- 
puscule isolé. Manuscrit reçu le 6 janvier 1959. 
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ÉTUDE DES CARACTÉRISTIQUES D'UN SPECTROGRAPHE £ 
UTILISANT UN CHAMP MAGNÉTIQUE A rH1 CONSTANT 


Par C. BASTARD, J. LAFOUCRIÈRE et R. MARGRITA, 


Institut de Physique Nucléaire de l’Université de Lyon. 


Résumé. — On étudie les caractéristiques d’un spectrographe 6 utilisant un champ magnétique 
inhomogène à rH constant. Une méthode d’approximation appliquée à l’ordinateur I. B. M. 650 
permet de construire les trajectoires gauches des particules chargées et de déterminer la forme des 


aberrations de sphéricité axiales. 


Abstract. — À study has been made of the characteristics of a B spectrograph using an inhomo- 
genous magnetic field of constant A. The twisted paths of the changed particles are computed 
with an I. B. M. 650 calculator for the determination of axial aberrations of sphericity. 


Nous avons vu, dans un précédent article [1], 
Pavantage que présentent les trajectoires tro- 
choïdales [2] pour l’études des spectres 8. Le champ 
magnétique utilisé est un champ de révolution 
possédant un plan de symétrie et tel que le pro- 
duit rH soit constant, r étant la distance à l’axe de 
révolution et H la valeur du champ dans le plan de 
symétrie. L'étude des trajectoires planes des parti- 
cules chargées, soumises à l’action de ce champ, nous 
a montré qu'un spectrographe $ utilisant ces tra- 
jectoires posséderait une focalisation radiale par- 
faite. Une équation approchée nous a permis de 
construire les trajectoires gauches, toutefois, nous 
avons dû nous placer dans les conditions de validité 
de cette équation, ce qui ne permettait pas d’uti- 
lser au maximum les possibilités du spectrographe. 
Nous avons repris l’étude des trajectoires gauches 
à partir des équations générales de manière à 
pouvoir discuter les caractéristiques du spectro- 
graphe dans le cas général. Signalons que les calculs 

-ont été effectuées à l’aide d’un ordinateur 
I. B. M. 650. 


I. Rappel des équations des trajectoires. —— Les 
particules de charge e décrivent des trajectoires 
dont les équations sont (1), (3) (coordonnées cylin- 
driques) 


mr” = mr p?—+ eg à (À) 


I { mr? o' + erA = C 
d 


LE 


mZ" = ep — (rA) 
où m est la masse relativiste et À la composante 
suivant Oz du potentiel vecteur, les deux autres 
composantes étant nulles. La fonction r A (rz) est 
égale à 
2\ 1/2 
Kr (1 se =. | (4) 


r 


avec À =r, H5, r, étant la distance de la source 
à Oz et À, la valeur du champ dans le plan de 
symétrie. 

La constante C peut s'exprimer de la manière 
suivante : appelons 0, l’angle de la vitesse ini- 
tiale v, avec un plan parallèle au plan de symétrie, 
ü, l’angle de la tangente au cercle de rayon rs, 
centré sur l’axe Oz avec la projection sur le plan 
de ce cercle de la tangente à la trajectoire, on a 


To Po = Lo COS 0, cos w 


d’où 
C = rifmus cos 0, cos v, + eK). 
Posons 
cK 
Kk— mt C = ro mvo (cos 0, cos us + k). 
0 £ 


Le système I s'écrit : 
lp" = volro (cos 0, cos u5 + k) — K(r? + z2)1/2] 
RATS kv, 
cle + #7) 


kvS oz 
74 Lt 0 
FR Eau 


Lorsque 0, — 0, les trajectoires sont situées dans 
le plan de symétrie et le mouvement est défini par 
les deux intégrales premières 


r'2 + r2 o’2 ae 32 — vè 
III} mr? o” + erA 2; Cr 


Intégrales qui s’intègrent facilement (3). 

Si # < 1 les trajectoires ont des branches infi- 
nies, Ce Cas ne nous intéresse pas. 

Si kZ1: on obtient des trajectoires tro- 
choïdales (fig. 1). 
. Quelque soit l'angle u,, toutes les trajectoires 
issues d'un même point S se recoupent sur le 


PR 
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cercle r 7, à une distance angulaire de $S égale 
à la précession 
DD TT (1 — : 


ea) 


IBtes 416 


On montre également que les trajectoires corres- 
pondant à 2 angles u, et u, se déduisent l’une de 
l’autre par une rotation et une homothétie de 
rapport 

cos u, + k 
COS Uo + k 


II. Caractéristiques du spectrographe. — Les 
caractéristiques du spectrographe vont dépendre 
de plusieurs facteurs à savoir 4, ro, Uo et 05, il faut 
tenir compte aussi des dimensions des trajectoires, 
ceci ayant surtout une importance en vue de la 
réalisation pratique. 


19 FOGALISATION RADIALE. — Nous avons vu 
dans le paragraphe précédent que les trajectoires 
situées dans le plan de symétrie et issues d’un 
point source S se recoupaient toutes en un même 
point situé sur le cercle r — r, à une distance angu- 
laire de S égale à 

k 
Ô —= DT (1 = el 
la focalisation radiale est donc parfaite et nous 
pouvons prendre une ouverture radiale importante. 
La direction moyenne de tir a été prise perpendi- 
culaire au cercle r = r,, cette disposition nous per- 
met en effet de séparer facilement la source et le 
dispositif récepteur par un écran (fig. 2a), celui-ci 
sert également à arrêter les trajectoires parasites 
qui correspondraient à des valeurs de Æ donnant 
des précessions sous-multiples de la précession, des 
trajectoires des particules détectées. Une direction 
moyenne de tir parallèle au cercle r — rs (fig. 26) 
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ne permettrait pas d’arrêter ces trajectoires. L’ou- 
verture radiale a été prise égale à 909, valeur qui 
nous a semblé la plus raisonnable pour placer 
l'écran précédent et le dispositif récepteur. 


HG? 


Deux directions de tir perpendiculaires au 
cercle r —7r, sont possibles, l’une dirigée vers 
l'intérieur du cercle r =.rj.(u, — 90°, fig 30), 
l’autre vers l'extérieur (4, — — 909 fig. 3b) nous 


1e 8 


verrons plus loin s’il existe un avantage théorique à 
choisir l’une plutôt que l’autre, nous pouvons 
cependant déjà remarquer que la première solution 
est la meilleure au point de vue pratique car elle 
laisse plus de place pour le dispositif récepteur. 


20 DIMENSIONS DES TRAJECTOIRES. — On peut 
montrer [3]: 
19 Que les trajectoires planes restent comprises 


entre deux cercles de rayons 


__ ro (COS uo + k) r 
PERS Le Cr 
ro (COS up + k) 


k 


Tm = 


étant le rayon moyen c’est-à-dire la valeur de r 
pour laquelle de/dr — 0. 

20 Queles trajectoires gauches restent situées à 
l'intérieur d’une surface de révolution d’axe Oz 
dont la méridienne a l’aspect de la figure 4. Les 
courbes tracées sont des ellipses dont le grand axe 


est égal à 
kr, (cos 0, cos us + k) 


k2—1 
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et le petit axe à 
ro (cos 00 cos uo + À) 


(x? RS 1)772 
Zi(cm) 
20 
10 
Ce (cm) 
û 40 
_10 
_20 
ETC 


Les parallèles dans le plan de symétrie ont pour 
rayons 


Re (cos 0, cos uo + À) 
2 D 0 


k LT 
et 
, ro (cos 6, cos us + k) 
RAS CE SAN RS 


Le tableau ci-dessous donne les valeurs de à, 
Ta le, Tm pour les valeurs de # allant de 1,3 à 2 et 
pour 4, égal à 450, 900 et 1359, r, ayant été pris 
égal à 20 cm. 


k Ô Ug — &59 Up — 900 Up — 1359 
ri = 17,5 cm r, = 11,3 cm r, — 5,2 cm 
198220902287 0133-8 cm7. 1060 cm7 50) 01cm 
Pre 0 ICE 20 C7 7-9 IRC 
rl 01C0 7 ETC ES SIC 
12002 7 T0" Siemer 0Rcme re 6m 
Tim — 00,1 C7» — 20m 0mern—/;9)ICmM 
r1 = 17,7 cm 7, — 12 cm r, — 6,3 cm 
IS PPS GE = CD MONTRE 07 on 
Fr 29) 0CMEr» 20cm, MD; 0ICM 
Pi 1/02 3 CMEr 6; 9m 
ONE TE CU ir er PIN CUT 
T0 28 8 0mern— 20 CNE: — 1412 1Cm 
ri 178 cmerr 12 EME EE CI 
LT RON En OO CN 7 = Po 
Tn 020) C7 220 Cl IC 
2 — 79cm 71-12 9 cer; 6 Cm 
PCR NM NO VAS on ZA = 2718 Gi 
Fm = 27,007 200Cm 77» —d2;hem 
F1 = 18 cm r,.— 13,1 cm 7, — 8,2 cm 
TM A TIC VA) din 7e = AE din 
Im 2 ACMer» — 20cm 2:01 
Tr, = 18 .cm r, — 13,3 cm r, — 8,6 cm 
25002007 = HE lCoEr-— 20 CmEr —#25 90cm 
Im 21; CNET 20cm M2 CM 
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dx 
est donné par la formule y Thu 


cisse curviligne du point image sur le cercle r =ro 


x étant l’abs- 


et oH égal à HET mv, étant l'impulsion initiale de la 
e 


particule 


les 


La figure 5 représente la courbe donnant la 
variation du pouvoir dispersif en fonction de k. 


49 FOGALISATION AXIALE. — L'étude de la foca- 
lisation axiale nous amène à construire les trajec- 
toires gauches solutions du système IT, ce système 
n'étant pas intégrable, nous devons utiliser une 
méthode d’approximation, les calculs ont été effec- 
tués à l’aide d’un ordinateur I. B. M. 650. La for- 
mule d’approximation utiliséeest une formule poly- 
nomiale due à Milne [5] elle est identique à celle 
déjà employée pour résoudre l'équation approchée 
des méridiennes (1). Cette formule fait intervenir 
une récurrence entre trois points successifs : 


Pn4ai = 2n — Fri + ra 

Zn+1 = 22n — Zn + 1274. 
Le premier point est obtenu à l’aide des formules 
n+i = 7h + Un (#+: = 2n + tan 


PTE Us 
PaEi= En À bn Horn Lines = 5 + + Que 
Une étude complète nécessiterait la construction 
des trajectoires correspondant à des valeurs de k 
allant de 1,3 à 2. 6, variant de 0 à 200 et w, de 450 


Ru IN° 8-9 


DAS étde 45021250 Cette étude conduisant 


à un calcul trop important, nous l’avons séparé 
en deux. La construction des méridiennes z — f(r) 
pour # — 2 à l’aide de l’équation approchée (1) met 
en évidence une concentration des trajectoires au 
voisinage du plan z = 0 mais pour une valeur de r 
différente de r,. Dans un premier calcul nous avons 


construit les courbes z —f(r) pour un angle 6, 


donné 9%, deux valeurs de u, : 900 et — 909, £ pre- 
nant les valeurs 1,3 ; 1,4 ; 1,5 ; 1,6 ; 2. Ceci afin de 
choisir la valeur de 4 et la direction de tir donnant 
la meilleure focalisation. Ce choix étant fait, nous 
avons construit les trajectoires correspondant à 
différentes valeurs de 0, et de us. 


III. Choix du paramètre k et de la direction de 
tir. — Les courbes de la figure 6 donnent la varia- 
tion en fonction de k de la cote z; du point d’inter- 
section des trajectoires avec le cylindre r =r, 
pour 49 —= 900 et uw, — — 900. Les deux courbes 
sont sensiblement identiques, il n’y a pas lieu de 
préférer l’une à l’autre au point de vue théorique, 
nous choisirons 4, — 909 car nous avons vu plus 
haut que cette valeur est plus satisfaisante que 
u — 909 du point de vue pratique. 


ICONS 


La courbe de la figure 6 correspondant à 
uo — 90° montre que pour la focalisation axiale 
deux valeurs de Æ sont favorables 4 — 1,25 et 
k — 2,04, théoriquement la première solution est 
la meilleure car elle donne un pouvoir dispersif 
bien plus grand que celui de la deuxième, par contre 
elle est pratiquement beaucoup moins satisfaisante 
par suite de l’encombrement des trajectoires 
comme le montre le tableau du deuxième para- 
graphe, le rayon maximum r, pour À —1,25 et 
uÿ = 45° est en effet égal à 156 cm, il serait difli- 
cile d'obtenir un champ en r-1 sur une aussi grande 
distance. D'autre part la cote maximum Z, des 
trajectoires est égale à 17,8 cm pour k — 1,25 et 
9,4 em pour k — 2,04, un entrefer de 36 cm ren- 
drait très difficile la production d’un champ assez 
intense. Pour # —1,25 la précession est égale 
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à 2400 ce qui oblige à utiliser une cuve à vide dont 
le volume est pratiquement égal à celui compris 
entre deux cylindres de rayons 160 em et 4 cm 
(rayons maximum et minimum des trajectoires 
planes) et de hauteur 36 cm soit environ 300 litres, 
par contre pour # — 2,04 on peut se contenter 
(voir fig. 1) d’une cuve dont le volume est celui 
compris entre deux cylindres de rayons 55 cm et 
8 cm, de hauteur 20 cm et deux plans diamétraux 
faisant un angle de 1209 soit un volume de 60 litres 
environ. | 

La réalisation d’un spectrographe avec £ = 1,25 
présenterait donc de grosses difficultés pratiques 
pour la production du champ magnétique et 
l’obtention du vide, nous avons donc effectué la 
2e partie du calcul avec 4 = 2,04 


IV. Étude détaillée dans le cas où k — 2,04. — 
1) FOGALISATION RADIALE. — Comme nous l’avons 
vu dans le paragraphe précédent la focalisation 
radiale est parfaite. La précession à est dans ce cas 
égale à 530. 


2) Pouvoir DisPERSIF. — Le pouvoir dispersif 
est donné par la formule (1) où 4 = 2,04et r, = 20 
on trouve 

AE 
Te 


Un spectrographe du type Siegbahn possé- 
dant une orbite stationnaire de 20 cm de rayon a 


Says 3 80 
un pouvoir dispersif égal à y =—,et un spectro+ 


o H 


. graphe à champ uniforme dont les trajectoires ont 


pour rayon 20 cm un pouvoir dispersif égal 


à Y DATE) il faut remarquer toutefois que dans 
notre cas les trajectoires s’écartent assez loin du 
cercle r —r, (voir fig. 1) et que l’encombrement 
est beaucoup plus important que celui d’un spec- 
trographe du type Siegbahn ou d’un spectrographe 


à champ uniforme. 


30 FOCALISATION AXIALE. — Nous avons cons- 
truit les trajectoires gauchesen utilisant la méthode 
de Milne à l'ordinateur 650, les valeurs de 0, et & 
utilisées sont les suivantes : 


uo — 900 0, — 3°, 60, 90, 120, 150, 180, 210 
Uo = 450 0, — 3°, 90, 150, 210 
Up = 135° 0, — 30, 90, 150, 240. 


Les courbes des figures 7 et 8 représentent la 
variation de z en fonction de r c’est-à-dire les 
courbes méridiennes des surfaces de révolution sur 
lesquelles sont tracées les courbes gauches (3) 
pour 0, égal à 30 et 9 et u4 égal à J00et 450. 

La variation de la cote z; du point d’intersection 
des trajectoires avec le cylindre r = rÿen fonction 
de 0, pour &, égal à 900, 450 et 1350 est représentée 
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sur la figure 9. En utilisant une équation appro- 
chée, nous avons trouvé pour #4 — 900 une fonc- 
tion linéaire en nous limitant aux valeurs de 0, 


to Z (cm) 
F8 


Es, 


inférieures à 109, il semble d’après la courbe précé- 
dente que la limite supérieure de 60, valable pour 
cette équation doit être prise un peu inférieure à 
cette valeur. 


Fic. 8. 


Pour u, — 1350 les trajectoires correspondant 
à 0, — 150et 0, — 210 ne recoupent pas le cylindre 
r — 20 em lors de la première boucle, ceci est dû au 


ZA (cm) 


| 
j 


USE 


5 


le 


3 


Fic. 9. 


fait qu’elles rencontrent les surfaces de révolution 
auxquelles elles sont tangentes (voir 1% para- 
graphe) à une distance inférieure à 20 em. La 
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“figure 10 représente la courbe z =/f(r) pour 


0, — 15 et la méridienne de la surface de révo- 
lution qui lui est tangente. 

Z'al(cm) 

10 

Peas r (cm) 
0 30 
10 
Fre-"0: 
49 FORMATION DE L'IMAGE D'UN POINT. — Un 


point source S donne une image sur le cylindre 
r — 20 cm. L’étude de la variation de la cote Z, 
des points d’intersection des trajectoires avec ce 
cylindre en fonction de la précession à permet de 
mettre en évidence la formation de cette image, la 
dimension de l’image dans la direction perpen- 
diculaire à Oz est en effet proportionnelle à à 
puisqu'elle est égale à r, à — 20 à, à étant exprimé 
en radian. 


ga 
20° 


15° 


0 1 L se 
SR 60 45° LORS 
Fi1G. 11. — Sur la coube médiane lire : U, = 90°. 


Lorsque l’angle O£est différent de 0 la précession 
n'est plus égale à 539 mais varie en fonction de 64. 


_ 


SUPREME 


1 


= 


4 


a — 
APTE 


SG SE 
; A (a 
dx 
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. La valeur dela précessionest donnée par l'intégrale 


de fe œ’(t) dt 
Tn 


où (1) est donné par le système II 


__ Lolro (cos 0, cos uo + k) — kr? + z2)1/2] 
Te 


(1) 


Cette intégrale a été calculée par la méthode de 
Simpson. Les courbes de la figure 11 donnent la 
variation de Ô en fonction de l’angle 6, pour 
Uo — 450, uy — 900 et u, — 1350. Les courbes 
de la figure 12 donnent la variation de z, en fonc- 
tion de à pour les mêmes valeurs de w,. 


12e (cm) Len OX MON PU SE Ed + 


rc 4? 


Si nous utilisons comme récepteur un compteur 
possédant une fenêtre de longueur L, nous voyons 
que pour une certaine valeur w; de 4, la courbe 
Z; — f(S) correspondant à cette valeur sera tan- 
gente à la droite z = L 

L'image d’un point source prend différentes 
formes suivant les valeurs de w,. 


Pour w, compris entre 45° et u; elle a l’espect de la 
figure 13a. 

Pour w, compris entre u; et 90° elle a l’aspect de la 
figure 130. 

Pour uw, compris entre 90° et 1359 elle a l’aspect de la 
figure 13c. 


En définitif, si w, varie de 900 à 1350 l’image a la 
forme de la surface hachurée de la figure 13d. 

Les résultats précédents montrent que la largeur 
de l’image correspond à une variation de la pré- 
cession égale environ à 2°; si on voulait détecter 
l’image en entier, il faudrait utiliser un compteur 
dont la fenêtre aurait 2,8 cm de long, il est évident 
que ceci entraînerait un mauvais pouvoir sépa- 
rateur, il est donc nécessaire de prendre-une largeur 
de fenêtre beaucoup moins importante, l’image 
d’un point se réduit alors pratiquement à une sur- 
face rectangulaire. 
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Dans un prochain article nous étudierons la posi- 
tion et la forme à donner aux diaphragmes, le 
pouvoir séparateur et la luminosité du spectro- 
graphe. 
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Nous tenons à remercier M. le Pr J. Thibaud, 
Directeur de l’Institut de Physique Nucléaire pour 
l'intérêt qu'il a bien voulu prendre à notre étude 
et les facilités qu'il nous a accordées. 

Manuscrit reçu le 11 mai 1959. 


BIBLIOGRAPHIE 


[1] BasrarD (C.) et LaroucriÈRE (J.), J. Physique Rad., 
1958, 19, 674. 

[2] Trisaup (J.), C. R. Acad. Sc., 1933, 197, 447, 915, 1629 
et Phys. Rev., 1934, 45, 787. 


[3] Basrarp (C.), Thèse, 3° cycle. 

[4] Moussa (A.) et LAFOUCRIÈRE ACER EACIAESC 
1951, 288, 239. 

5] Mae (We), Numerical calculus 142. 


LE JOURNAL DE PHYSIQUE ET LE RADIUM TOME 20, AOUT-SEPTEMBRE 1959, PAGE 742. 


VYARIATIONS DES PROPRIÉTÉS OPTIQUES DE COUCHES MINCES D'OR 
EN FONCTION DE LEUR VITESSE DE FORMATION 


Par R:<PHILIE: 


Laboratoire de Physique Générale, Faculté des Sciences de Marseille. 


Résumé. — On a étudié l'influence, sur les propriétés optiques de couches minces d’or obtenues 
par évaporation, de la vitesse de formation du dépôt. On a utilisé pour cela trois séries de couches 
préparées respectivement avec des vitesses de formation de 0,8, 3 et 25 my à la minute. Toutes 
les couches d’une même série, dont les épaisseurs variaient entre 0 et 70 mu, ont été préparées au 
cours de deux évaporations conduites pratiquement à la même vitesse. On montre que les pro- 
priétés optiques des couches varient beaucoup avec la vitesse de projection. En particulier l’absorp- 
tion À’ d’une couche d’or de 60 my d'épaisseur peut être multipliée par 25 lorsque la vitesse de 
formation passe de 25 à 0,8 mu à la minute. 


Abstract. — A study has been made of the influence of the speed of formation on the optical 
properties of evaporated gold films. Three series of evaporated layers have been prepared res- 
pectively with a speed of 08, 3 and 25 my by minute. Evaporated films of similar series, the 
thickness of which varied from 0 to 70 mu, have been prepared during two evaporations made as 
nearly as possible at the same speed. It was found that the optical properties vary a great 
deal according to the speed of formation. Specially, the absorption À’ of à 60 mu thick film 
may be multiplied by 25 when the speed of formation is increased from 25 to 0.8 mu by minute. 


I. Introduction. — Dans un article publié en 
1950, Sennet et Scott [1] ont étudié qualitati- 
vement au microscope électronique l'influence de la 
vitesse de formation sur la structure de couches 
minces métalliques obtenues par évaporation ther- 
mique. Ils ont ainsi montré que les couches métal- 
liques minces étaient d'autant plus lacunaires que 
leur vitesse de formation était plus faible. Ils ont 
ensuite étudié quantitativement les variations des 
propriétés optiques de ces couches minces avec leur 
structure. Mais leurs mesures sont incomplètes (ils 
ne mesurent pas le facteur de réflexion R’ à l’inté- 
rieur du support), peu précises ( 2 %) et il n’est 
tenu compte de la présence du support de la couche 
mince que par une correction «approximative ». Ces 
mesures n’ont donc pas pu conduire à la déter- 
mination des constantes optiques des couches. 


IT. Mesures. — Dans le présent travail nous étu- 
dions 3 séries de couches minces d’or déposées sur 
un support en quartz, sous un vide voisin de 
1075 mm de mercure, avec des vitesses de for- 
mation de 0,8, 3 et 25 mu à la minute. On appelle 
vitesse de formation d’une couche le quotient de 
son épaisseur « massique » par la durée de l’évapo- 
ration (1). 

À l’aide d’un disque, genre disque de Talbot, nous 
avons pu obtenir simultanément, au cours de cha- 


(*) L’épaisseur « massique » est obtenue en pesant, avant 
et après l’évaporation thermique une lamelle couvre-objet 
de microscope de surface connue et en supposant que la 
densité du métal en couche mince est la même que celle du 
métal massif. On sait que, cette hypothèse est inexacte et ne 
constitue qu’une approximation commode. 


cune des évaporations, 9 couches minces d’épais- 
seurs croissant en progression arithmétique et 
formées rigoureusement dans les mêmes conditions 
(en particulier même vitesse de formation). Cha- 
cune des courbes présentées ci-dessous a donc été. 
obtenue à l’aide de deux évaporations seulement 
conduites à des vitesses de formation sensible- 
ment égales. 

Pour chacune des couches minces étudiées nous 
avons mesuré, à moins de 0,5 % près, le facteur de 
transmission 7, les facteurs de réflexion côté air À 
et côté support À”; nous en avons déduit les absorp- 
tions correspondantes À = I — R — T et 
A" —I— R'— T. L’évaporateur [2] ainsi. que le 
principe des mesures [3] ont déjà été décrits. Les 
courbes ci-dessous résument les résultats obtenus. 
L’auteur tient à la disposition des lecteurs qui 
voudraient en prendre connaissance les tableaux 
numériques des résultats expérimentaux obtenus. 


IIT. Discussions. — A) Les courbes des figures 1 
à 9 présentent, pour les 3 vitesses de formation 
étudiées, les variations de T, R, R', A et A! en 
fonction de l’épaisseur de la couche métallique. 


19 COURBES RELATIVES AU FACTEUR DE TRANS- 
MISSION T. — Pour toutes les longueurs d’onde 
étudiées et pour une épaisseur quelconque com- 
prise entre 0 et 70 mu, le facteur de transmission T 
est toujours d'autant plus grand que la vitesse de 
formation des couches est plus grande. C’est ainsi 
(fig. 1) que pour la longueur d’onde 6 438 À et 
pour une même épaisseur de 16 mu, T passe de 0,18 
pour une couche formée à la vitesse de 0,8 mu à la 
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Fic. 1,2 et 3. — Variations du facteur de transmission T 
et du facteur de réflexion côté air À, en fonction de 
l'épaisseur massique de la couche mince d’or, pour les 
3 vitesses de formation étudiées et pour les 3 longueurs 
d’onde 6 438 À ; 2 967 À ; 3 341 À. 
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minute à 0,51 pour une couche formée à la vitesse 
de 25 my à la minute. 

Cet écart est plus petit pour les radiations du 
spectre ultraviolet (fig. 2) que pour celles du spectre 
visible. On remarquera aussi que la courbe relative 
à T présente, surtout pour les plus fortes longueurs 
d'onde, un palier plus ou moins marqué suivant la 
vitesse de formation des couches minces. Ce palier 
est très marqué pour la radiation de longueur 
d'onde 6 438 À et pour une couche formée à la 
vitesse de 25 my à la minute ; T conserve alors la 
même valeur pour des couches minces d’épaisseurs 
comprises entre 5 et 12 mu. Pour les couches 
formées aux vitesses de 0,8 et 3 my à la minute, 
un léger palier peut encore être observé pour T mais 
uniquement pour des épaisseurs comprises entre 4 
et 5 my ; aussi passe-t-1l assez souvent inaperçu. 


20 COURBES RELATIVES AUX FACTEURS DE RÉ- 
FLEXION À ET R'. — Pour toutes les radiations 
étudiées ici le facteur de réflexion côté air R est 
d’abord, pour les faibles épaisseurs, d'autant plus 
grand que la vitesse de formation des couches est 
plus petite, la variation de À avec cette vitesse 
reste alors assez faible (0,03 à 0,05). Au delà d’une 
certaine épaisseur, variable avec la longueur d’onde 
(24 mu pour à — 6 438 À et À — 4 046 À ; 10 mu 
pour À — 3 341 À et À — 2 967 À) c’est l'inverse 
qui se produit. La couche métallique d’or réfléchit 
d'autant plus qu’elle a été plus rapidement formée ; 
la variation de À avec la vitesse d’évaporation 
devient alors très sensible. On trouve, par exemple, 
pour une même épaisseur de 60 mu et pour la lon- 
gueur d'onde 6 438 À, que R passe de 0,6 à 0,85 
lorsque la vitesse passe de 0,8 à 25 my à la minute 
(fig. 1). Pour la même épaisseur et les mêmes 
vitesses d’évaporation À passe de 0,20 à 0,35 pour 
la radiation de longueur d’onde 2 967 À (fig. 2). 
On remarquera aussi que, pour À — 2 967 À, R pré- 
sente un léger minimum pour une épaisseur voisine 
de 10 muet pour la plus faible vitesse de formation. 

Pour les autres radiations on observe soit un 
palier (pour À —3 341 À, fig. 3), soit une simple 
variation de pente, surtout marquée pour les radia- 
tions ultraviolettes et pour les vitesses de for- 
mation très lentes. 

Le minimum du facteur de réflexion côté sup- 
port R’ ainsi que l’épaisseur correspondante de la 
couche mince varient avec la vitesse de formation 
(fig. 4). Ces variations du minimum de À’ restent 
toujours inférieures à 0,005, c’est dire qu'elles 


Nota. — Pour toutes les courbes ci-jointes nous adoptons 
les conventions suivantes : 


pour les couches d’or formées à la vitesse 
de 0,8 my à la minute ; 
—,—.—.— pour les couches d’or formées à la vitesse 
de 3 my à la minute ; , 
pour les couches d’or formées à la vitesse 
de 25 mu à la minute. 
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Fic. 4, 5,6 et 7. — Variations du facteur de réflexion côté 
support À’ et du coefficient d’absorption 4’, en fonction 
de l’épaisseur massique de la couche mince d’or, pour 
les 3 vitesses de formation étudiées et pour les 4 lon- 
gueurs d'onde ; 2 967 À ; 6438 À ; 5 461 À ; 3 341 À. 


nécessitent pour être mises en évidence, des mesures 
très précises. Là encore, pour les faibles épaisseurs 
(inférieures à 10 ou 20 my suivant la longueur 
d'onde), ce sont les couches formées le plus rapi- 
dement qui réfléchissent le moins, par contre, pour 
des épaisseurs supérieures c’est l’inverse qui se pro- 
duit. La variation de À’ peut alors devenir très 
grande ; c’est ainsi que pour une couche mince d’or 
de 60 my d'épaisseur, et pour À — 6 438 À, R' passe 
de 0,31 à 0,90 lorsque la vitesse d’évaporation 
passe de 0,8 à 25 my à la minute (fig. 5 et 6). 

On remarquera enfin que les facteurs de 
réflexion À et R' peuvent présenter un maximum 
d'autant plus net que les couches sont formées plus 
lentement. On constate en particulier (fig. 7) pour 
la radiation de longueur d’onde 3 341 À, et pour 
une vitesse d’évaporation de 0,8 my à la minute, 
qu'une couche d’or de 60 my d'épaisseur réfléchit, 
côté support, la moitié moins (R' — 0,07) qu’une 
couche d’or de 30 my d'épaisseur (R! — 0,15). 


30 COURBES RELATIVES AUX ABSORPTIONS À 
ET À’, — On constate, pour une épaisseur donnée 
et pour toutes les longueurs d’onde étudiées, que 
les absorptions À et 4’ sont d'autant plus grandes 
que la vitesse de formation de la couche est plus 
petite. C’est ainsi que pour une épaisseur de 60 my 
et la radiation 6 438 À, A passe de 0,07 à 0,39 si la 
vitesse de formation de la couche passe de 25 my 
à 0,8 my à la minute (fig. 8). L’absorption A’ varie 
encore davantage puisque dans les mêmes condi- 
tions elle passe de 0,025 à 0,67 (fig. 5). Cet écart 
diminue avec la longueur d’onde de la radiation 
utilisée ; pour la même épaisseur et les mêmes 
vitesses de formation que précédemment mais pour 
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Fic. 8 et 9. — Variations du coeflicient d’absorption À, 
en fonction de l’épaisseur massique de la couche mince 
d’or, pour les 3 vitesses de formations étudiées et pour 
les 2 longueurs d’onde 6 438 À ; 2 967 À. 


la longueur d’onde 2 967 À, A passe de 0,60 à 0,78 
(fig. 9) tandis que A’ passe de 0,68 à 0,91. 

On constate de plus que les courbes relatives aux 
absorptions À et À’ ont, pour toutes les longueurs 
d'onde inférieures à 5 000 À, sensiblement la même 
allure générale (fig. 7, 9) et ne présentent pas de 
maximum, quelle que soit la vitesse de formation 
des couches. Par contre, pour les radiations de lon- 
gueurs d'onde supérieures à 5 000 À, ces courbes 
peuvent avoir des allures totalement différentes ; 
elles peuvent en effet, suivant la vitesse de for- 
mation des couches, présenter ou non un maximum 
plus ou moins accentué (fig. 5 et 8). 

B) Les courbes des figures 10 à 14 montrent les 
variations de 7, À, R', A et A’ en fonction de la 
longueur d’onde, pour les 3 vitesses de formation 
étudiées et pour des épaisseurs de 2 ; 10 ou 60 my. 
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On voit très bien sur ces courbes que, pour une 
épaisseur donnée, les couches formées à vitesse 
lente sont toujours les plus absorbantes (fig. 13 
et 14 tandis que les couches évaporées le plus 
rapidement ont toujours le plus fort facteur de 
transmission (fig. 10). Par contre, pour les facteurs 
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et pour des couches d’or d’épaisseurs massiques égales 
à 2 ; 10 ou 60 mu. 
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de réflexion À et R' on constate sur ces courbes 
que pour une épaisseur de 10 my ce sont les couches 
formées lentement qui réfléchissent le plus, tandis 
que pour une épaisseur de 60 my les couches 
formées lentement sont celles qui réfléchissent le : 
moins (fig. 11 et 12). 

Les couches minces d’or d'épaisseur 10 mu pré- 
sentent un maximum pour T (vers 5 200-5 300 À) 
et un minimum pour À et R' (vers 5 000 À) ainsi 
que pour À et 4’ vers 5 200-5 300 À. Pour des 
couches d’or de 60 my d'épaisseur seul subsiste le 
maximum de 7”. La position de ces extrêma varie 
du reste avec l’épaisseur de la couche considérée. 

La relation de Wolter (n A — À' = 0) elle-même 
subit l'influence de la vitesse de formation des 
couches mais surtout pour les très faibles vitesses. 
On sait que cette relation est moins bien vérifiée 
pour les radiations du spectre ultraviolet que pour 
celles du spectre visible. Nos résultats montrent 
que, dans le spectre visible, pour des couches d’or 
formées à la vitesse de 0,8 mu à la minute, cette 
relation n’est vérifiée que pour des couches d’épais- 
seurs inférieures à 5 ou 6 mu. Par contre, pour des 
couches d’or formées à la vitesse de 3 ou de 25 mu à 
la minute cette relation est vérifiée dans le spectre 
visible pour des couches d’épaisseurs plus grandes. 


Conclusion. — Les propriétés optiques d’une 
couche métallique d’or d'épaisseur donnée peuvent 
être très différentes suivant la rapidité avec laquelle 
la couche a été formée. Bien d’autres facteurs 
peuvent aussi intervenir, ce qui explique le manque 
de concordance entre les résultats des divers 
auteurs et la nécessité pour chaque expérimen- 
tateur de préciser les conditions dans lesquelles il a 
opéré. 


Manuscrit reçu le 10 avril 1959. 
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UNE VERSION MODIFIÉE DE LA THÉORIE DU CHAMP UNIFIÉ D’EINSTEIN 


Par Jack LÉVY, 


Institut Henri-Poincaré, Paris. 


Résumé. — Pour éviter les difficultés auxquelles se heurtent les théories usuelles du champ 
unifié, nous sommes conduits à introduire dans un lagrangien un terme quadratique par rapport 
au vecteur de torsion de la connexion affine. Cette introduction est justifiée par les résultats 
de la théorie de M. H. Eyraud et un schéma hydrodynamique de la théorie des électrons de Dirac. 


Abstract — To avoid restrictions and difficulties encountered in classical theories of Einstein’s 
unified field, we are lead to introduce in a Lagrangian formalism, a tensor built as T4 , from the 


torsion vector. 


Fields equations in two equivalent sets are deduced, and a current is pointed 


out, able to give non-vanishing equations of motion. 


Introduction. — L’unification des champs gra- 
vitationnel et électromagnétique peut se concevoir 
dans le cadre d’une géométrie affine quadridi- 
mensionnelle. Sur un réseau curviligne x" on attache 
des repères locaux e, en des points infiniment 
voisins. La connexion de deux tels systèmes de 
repères fournira le type de géométrisation. Une 
extension directe de la relativité générale donne la 
formule de liaison 


des == Le dx? ey y 


où les [°,> sont les coefficients de connexion affine, 
fonctions à propriétés de symétrie non particu- 
larisées. 

À partir de cette distribution fondamentale non 
tensorielle, on peut par des intégrations locales [1] 
construire des champs tensoriels, courbures, tor- 
sion. Muni de cette surstructure, on peut alors 
former des lagrangiens, en déduire des équations 
de champs covariantes et compatibles. 

Après un bref rappel des théories usuelles et de 
leurs difficultés, nous proposerons un moyen de 
remédier à ces dernières. Nous justifierons et pré- 
ciserons le choix de notre lagrangien par un rap- 
prochement avec deux théories, l’une de M. Eyraud, 
Pautre de Dirac. 


A. Théories usuelles. — 1. LA THÉORIE MIXTE 
d’Einstein juxtapose un nouveau champ tensoriel 
défini par le tenseur métrique g*’, aux tenseurs 
directement construits en géométrie affine. On 


posera :. 


qu = Viglgu; gdet. de guy;  Ruv tenseur de Ricci: 


La comparaison avec la relativité générale 
suggère alors pour un lagrangien le choix L, ou Lo, 


dont on déduit par variations libres et indépen- 
dantes des g,, et [8 les équations de champ J,. 


Lo = qu" Roy JE Æ au Rav — 27 VTgl 


| (a) 9 3 p = de gé— ge Dés + géo Déc + 9° Dé 


2 v 
Lo M 0) 
(b) 29 = 0 
(c) 


Ruv À Suv- 


À Ne EF) 0 


introduite par Einstein définit un système « fort ». 
Schrôdinger, par l’intermédiaire du changement de 


La condition arbitraire PF, = 


2 . à 
connexion Le —17,$ +3 6 l,, déduit un système 


«faible » d’un principe variationnel. 

Quoi qu'il en soit, le vecteur de torsion est inac- 
cessible, l'identification des divers tenseurs malaisée 
et les identités de conservation ne permettent pas 
de retrouver les lois du mouvement d’une particule 
chargée dans un champ électromagnétique. Gette 
impossibilité est sans doute due ([5] p. 120 et [61) 
à l'équation Z,(b) par là trop restrictive. 


2. Des compLicaTions du tenseur de Ricci ([5] 
p. 129 et [7]) permettent de s’en affranchir, mais 
il semble plus simple, à partir d’une connexion 
affine très générale, d’évaluer systématiquement les 
ressources brutes de notre géométrie en tenseurs 
du second rang. 


3. PAR coNTRACTION du tenseur de Riemann 


RévolT) = do TÉv — ds ne + 1e Te == ne I 
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on obtient : 
Ron Lie | PS en (3) 
La torsion fournit les tenseurs : 
Tu IV du Ly — dv Ty. (3’) 
En introduisant la notion de connexion trans- 


posée D, — lé, on obtient RE) et Pa) 

On peut alors former une combinaison linéaire 
de ces divers tenseurs et chercher les équations de 
champs pour le lagrangien : 


L = ge[Ru(T) + aRyT) SF bPy(T) DS cPyy(T) 
M 'aCT 
On trouve en particulier [8] une généralisation 
de Z,b 
A De ge = Ba ge TÉs + [20 ge + : di + a) g®| T, (4) 


A = 4b + ce 1 


avec a(b + 4c — a). 


On voit alors que sans altérer beaucoup le prin- 
cipe de cette généralisation, on peut en simplifier 
la forme en supposant b et c quelconques, a en d 
nuls. Le cas ou d £ O semble plus favorable puisque 
ainsi la notion supplémentaire de transposition 
n’est pas introduite, et une analogie apparaît avec 
deux théories : l’une est due à M. H. Eyraud [9] 
et par la simple introduction du terme l, F, dans 
un hamiltonien du type : 


A — Vo Ron 2 B(du Ty, — d, Tu) Se Yu P,| 


permet la définition d’un courant ; l’autre théorie, 
euclidienne, reconsidère le formalisme électro- 
magnétique classique et aboutit à un schéma 
hydrodynamique relativiste. 


4. LA NOUVELLE THÉORIE CLASSIQUE des élec- 
trons de Dirac introduit en effet une condition non 
linéaire sur les potentiels en remplaçant la conser- 
vation 2, À — 0 par la normalisation À, A —#?. 
Partant du lagrangien 


Re — 7 Fay Fuv + : MA Au — k?) 
on peut alors définir une proportionnalité courant 
potentiel 

de Fue — À AU, 

Dirac montre qu’un mouvement est possible sui- 
vant les équations de Lorentz (4 — masse/charge). 
En effet, soit À, le potentiel en absence de matière 
chargée (À — 0), par la création de charges petites, 
le potentiel varie peu et reste normalisé : 
(Au + du 5) (4 + 25) = ne. 
L’impulsion est égale à : 


Ay = x + du S 


I, = eds. (5) 
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Dans une seconde étape Dirac considère un flux 
d'électricité à champ continu de vitesse v,. 


et remarque que la loi de force de Lorentz peut 
s’écrire successivement : 


dos 


k AL Ue Feu. | (7) 
IIue pv, = 0 (8) 
avec 
Lluv = Fuv + (y Vu — du Uy) = du 1, — dv lu (9) 
1, = Ag — Ky. (9°) 


Ces grandeurs peuvent être considérées respec- 
tivement comme tenseur tourbillon et vecteur 
impulsion d’un fluide parfait chargé [14, 15]. 

Soit alors l’espace vectoriel x à deux dimensions 


des vecteurs V (tes les vecteurs-vitesse v, et tour- 


billon + — ; ne 4 ra) pour lesquels : 


IL 149 —10; 


Ces vecteurs dits de quasi-tourbillon, peuvent 
être orientés dans le temps, l’espace ou être iso- 
tropes. Soit encore un système de couples de 
fonctions E* et n, (a —1, 2 ...) constantes sur 7. 
On a 


US do E% — T9 do NE —= ... == (]}, (10) 
Donc d’après (9) on peut écrire le tenseur tour- 


billon sous la forme 


Iuy = (Ou E% dv na — dy E€ du Na) 
et l'impulsion 


Lu = Ay — kvy = Ea dy Na (11) 
ce qui correspond au mouvement tourbillonnaire 
du flux d'électricité. 

Ces équations sont déductibles d’un lagrangien 
où Ë, », Vy, A SOnt variés : 


1 il 


à Ne F,, étant fonction de £, n d’après (9) et 

Ces résultats nous suggèrent l'introduction expli- 
cite en théorie unitaire d’une condition de norma- 
lisation, en espérant par là, fournir une interpré- 
tation physique des potentiels, un schéma hydro- 
dynamique et en précisant ce cadre, la signification 
des divers tenseurs. 


\ 


B. La condition de normalisation en théorie 
unitaire. — 1. ÉQuATIONS DES cHAMPs. — Nous 
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adopterons l'expression suivante du lagrangien (1). 
\ 

; © = gUVERT) + mÜSu,s — Syu)] 

+ QT, — Au) + a(gu Su 5, — 2 VTg] a). 


En effet, en jouant sur leterme A#(T, — k,) nous 
pourrons faire porter la condition de normalisation 
sur un vecteur quelconque S, et non plus néces- 
sairement sur le vecteur de torsion [L', ; ceci afin 
d'éviter toute identification ou interprétation anti- 
cipée des champs. Le terme en a impose la norma- 
lisation. Nous verrons par la suite la nécessité du 

terme en m ; R,, provient de l'extension directe de 
la relativité générale. Précisons les notations : 

R,(l) tenseur de Ricci de la connexion quel- 
conque [”,4 ; 

W,,(L) qui est, lui, construit sur la connexion L,£à 
vecteur de torsion nul ; 

A —\/]g] AY, a scalaire, multiplicateurs de 
Lagrange ; 

k, vecteur susceptible d’être fonction de S,,. 

Les variations portent sur g#, [4 S,, X4, a, et 
donnent les systèmes d'équations : 


BY 
I CRE L— gÊY 1 —° + ge les 
L de gPe + guy LÉ, _. : gfe I, + AB — 0 


I, d ge = Sur 


IT Ru + mSus — Syu) + aSu Sy = à à? guv 


III age Se — do ge? (5 DE + m 5) 


IV Ro Ty 
V gum Su 5 — 2 V'Ig] «?. 


On voit alors le caractère nécessaire de l’en- 
semble des champs auxiliaires A, a, m, «. Le vec- 
teur A# par I, dote la théorie d’un «courant » dont 
nous connaissons la nécessité dans la détermination 
du mouvement des particules. D’où l’exigence 
d’après la signification de Y, dans (1) d’une condi- 
tion sur le vecteur de torsion, puisque l’existence 
de Y, entraîne l’équation IV. Cette condition est 
ici explicitée par a et k, en une normalisation. 
L'examen des équations I à V montre d’autre part 
que : 

— %% étant nécessairement différent de zéro 
pour s'affranchir de 1,4. 

— Sia = 0,IIlet I, conduisent à 


es [2 2 ke go) _ 
do 8 (536) + 20m 8€) F 


LE + m D do QE + © do m = 0 
œ 


(4) Pour ce choix et le détail des calculs voir [8]. 
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c’est-à-dire à I, avec les hypothèses que nous 
serons rapidement amenés à faire sur k, et m. 

— Sim — 0 on doit avoir 4, £ 0 ce qui élimine 
de la théorie la possibilité de retrouver un système 
faible à vecteur de torsion nul. 

Nous allons maintenant supposer successi- 
vement #, nul et k, identique à kS, Il vient les 
deux systèmes. 


d(7 g%) —= À g— So 
? 
de QËe + Qu" Lis + 3 de QE = 0 


II Wuv(L) +- MS uv ES Sy, 1) + aSy Sye— au? Suv 
IV Ps) | 


(a QUY Su Sy = 202 V[g]. 


Ce système À sera dit faible (Z, — 0). Le vecteur 
de torsion a disparu, la géométrie n’est plus géné- 
rale, mais on a dû introduire pour remédier à cette 
restriction, un vecteur $,. Introduction phénomé- 
nologique à rapprocher de celle du tenseur g,,. 

Le système B homogène voit l'élimination de A# 
et 5, et s’écrit dans la connexion générale l'£, : 


BY 2 2 alk 
I res — Y ADN (e # 
Perso ler Sc 
ÉgeT, 
I ap alk ao 
NN cire, LE 


7) 
T5 De gÉe + ge TN + à ge To +50 gŸ — 
II Ryv(T) + mk(Lu,v = Tu) + ak? Te ES — à à? &uv 
IIL guy Tu Ly = 202/k2.V [gl]. 


Remarquons les valeurs interdites pour m: 


2k 
Ma FO et M Re 
2. ÉQUIVALENCE DES SysTÈMES. — Le sys- 


tème B absorbe le champ S, dans le schéma de 
connexion affine qui supporte alors la condition de 
normalisation. C’est donc un système plus unifié 
que À qui comporte en plus un vecteur S, ; mais en 
apparence seulement car nous allons montrer que 
ces deux systèmes sont transformés l’un de l’autre 
par une opération sur la connexion affine, opération 
qui n’affecte pas l’équation des géodésiques. 

Soit en effet un À-transformation définie [3, 4] 
par un vecteur À, : 


PT SN. 


Elle laisse bien invariantes [2, 3, 4] les équations 
géodésiques telles qu’Eisenhart les a données HOT: 
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On peut d’ailleurs rapprocher ces transformations 
des déformations introduites par M. Eyraud [9]. 

Soit alors 5%,, $, une connexion et son vecteur de 
torsion. Posons : 


Lér = S6y + 36 A7 + DY Be. 


Portons dans l’équation I du système À, on 
retrouve I de B, à condition de poser [8] : 


2 
Aa EN 152 


2k 
5 MA U— k (rer +7) 


Bei—00 
On retrouve II, avec 


9 
mg k = M4 — =: 
: 3 


Ce qui nous conduit à prendre 4? = 1 et 
MA + MB = 3 


Les systèmes À et B sont donc mathémati- 
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quement et physiquement équivalents, en ce sens 
qu'une transformation conservant le parallélisme 
échange les équations de champs déduites de 


Sa = QUIWuv(T) + Mma(Suy — Svu)l 
+ QU Lu + alge Su Ss — 2 V [g] a?) 


et celles données par : 


Lp = QUIRa(T) + mp Tu,v — Lvu)] 


+ a(guT, T, — 2 V/{g] «?). 


L'un des systèmes, À, est dégénéré de l’autre B 
par un choix de « jauge » convenable. 

Ainsi nous avons levé une restriction de la théorie 
du champ unifié en douant l’équation en diver- 
gence d’un second membre non nul par l’intro- 
duction dans le lagrangien du tenseur F, F, sous | 
forme d’une normalisation. 


Manuscrit reçu le 5 mai 1959. 
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ÉTUDE'DES TRANSITIONS ENTRE ÉTATS MÉTALLIQUES ET ISOLANTS. 
APPLICATION AUX BANDES D'IMPURETÉ ET AUX ANTIFERROMAGNÉTIQUES. II. 


Par J. Des CLOIZEAUX, 


École Normale Supérieure, Laboratoire de Physique. 


Résumé. — Dans certaines conditions un gaz d'électrons en mouvement dans un solide peut 
présenter un caractère isolant en contradiction avec les prévisions de la théorie des bandes. Un 
modèle simple a été proposé précédemment. Une étude plus précise corrobore les résultats obtenus 
et montre qu’à côté des excitations individuelles, il existe deux types d’excitations collectives, les 
plasmons et les spinons. L’application du modèle aux systèmes désordonnés est particulièrement 


discutée. 


.… Abstract. — Under definite assumptions, an electron gas in a solid may behave as an insulator, 
in contradiction to the predictions of band theory. À simple model has been proposed pre- 
viously. A more careful study gives the same results but shows that besides individual exci- 
tations there are two kinds of collective excitations : plasmons and spinons. An application of the 


model to disordered systems is discussed. 


VIT 


EXCITATIONS INDIVIDUELLES 
ET EXCITATIONS COLLECTIVES. 
THÉORIE PLUS COMPLÈTE. 


1. Généralités. — Jusqu'ici [28] nous avons 
négligé délibérément dans l’hamiltonien H, tous 
les termes non diagonaux par rapport aux états A%, 

% A et Bf Examinons mintenant les modifi- 
cations qui sont introduites dans notre modèle, 
quand on veut tenir compte de ces termes. L'état 
de base alors n’est plus F, mais un état légèrement 
modifié F;. 

Or, au paragraphe III-7, nous avons dressé la 
liste des opérateurs S,4 et S,, associés aux exci- 
tations possibles de l’état de base F, de l’hamil- 
tonien diagonal. Il est alors facile de voir que, si 
l’on veut pousser l’approximation d’un cran plus 
loin, il faudra extraire de l’hamiltonien non dia- 
gonal y» les termes linéaires ou bilinéaires par 
rapport aux opérateurs Si et S,x ; Soit Hxypy cette 
fraction utile de l’hamiltonien. Le problème con- 
siste alors à calculer les énergies et les excitations 
d’un système dont l’hamiltonien est maintenant : 


Ho = Hp + Hnpu. 


Ce problème a été résolu par Sawada et Sawada, 
Brueckner, Fukuda et Brout [4] pour un gaz 
d'électrons libres au zéro absolu (en ne tenant 
compte que des excitations qui conservent le spin, 
et sont seules importantes aux fortes densités). 

La méthode utilisée par ces auteurs permet à la 
fois d'évaluer l’énergie de corrélation de l’état de 
base et d’étudier la nature des excitations élémen- 


taires du système. En particulier, ce traitement 
rend compte de l’existence d’excitations de carac- 
tère collectif ou individuel, et permet de calculer 
les énergies d’excitation. 

Notre problème peut être résolu de façon ana- 
logue, bien qu’il soit un peu plus général. En effet 
à côté des excitations qui conservent le spin et qui 
conduisent à des excitations collectives du type 
plasmon, nous introduisons des excitations qui ren- 
versent le spin et conduisent, comme nous le 
verrons, à des excitations collectives du type 
spinon. Par ailleurs, pour étudier les diverses exci- 
tations du système, nous nous placerons à la tempé- 
rature T, cette hypothèse n’introduisant pas de 
difficultés supplémentaires et ne changeant rien au 
formalisme utilisé. 

Le but de ce calcul n’est pas d’évaluer l’énergie 
de corrélation. En effet, dans notre modèle, il y a 
séparation plus ou moins effective des électrons de 
spins opposés sur deux sous-réseaux ; l'énergie de 
corrélation a donc déjà été en grande partie absor- 
bée dans les termes qui proviennent de l’hamil- 
tonien diagonal. Comme nous le verrons, tout se 
passe comme si notre transformation antiferro- 
magnétique avait écranté le potentiel répulsif entre 
électrons. Par ailleurs, l'énergie de corrélation peut 
être évaluée par un simple calcul de perturbation 
au deuxième ordre, l'existence d’une énergie d’acti- 
vation e, assurant de toutes façons la convergence 
de ce calcul. Remarquons que, pour un gaz d’élec- 
trons libres, il n’en est pas de même, le calcul de 
perturbation du deuxième ordre diverge, et c’est 
justement cette divergence qui oblige à avoir 
recours à des procédés plus raffinés. : 

En ce qui concerne notre modèle, l'intérêt de la 
méthode réside plutôt dans l'étude des excitations. 
Il est, en effet, important de montrer que l'intro- 
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duction de termes non diagonaux dans notre hamil- 
tonien, ne détruit pas le caractère isolant du sys- 
tème. Les conclusions de cette étude sont alors les 
suivantes : 


a) Toutes les excitations correspondant au vec- 
teur d’onde gq, sont divisées en deux groupes dis- 
tincts, les excitations longitudinales qui conservent 
les spins et les excitations transversales qui les 
renversent. 


b) Les excitations longitudinales associées au 
vecteur d'onde q comportent un grand nombre 
d’excitations individuelles (5 — 
tation collective ou plasmon. Les énergies asso- 
ciées aux excitations individuelles de l’hamiltonien 
non diagonal sont sensiblement les mêmes que les 
énergies d’excitation correspondantes de l’hamil- 
tonien diagonal (le déplacement énergétique est de 
l’ordre 1/N). 

Le plasmon, quand il n’est. pas absorbé dans le 
fond continu des excitations individuelles, possède 
une énergie nettement supérieure à celle de toutes 
les excitations individuelles. 

c) Les excitations transversales de vecteur 
d'onde q comportent de même un grand nombre 


1) et une exci- 


2 
tation collective ou spinon. Les énergies associées 
aux excitations individuelles de l’hamiltonien non 
diagonal sont encore sensiblement les mêmes que 
les énergies d’excitation correspondantes de l’ha- 
miltonien diagonal (le déplacement énergétique est 
de l’ordre 1/N). 

Le spinon, quand il existe, possède une énergie 
inférieure à celle de toutes les excitations indi- 
viduelles. Si [q| —> 0, l'énergie du spinon tend vers 
Zéro. 

En résumé, pour toutes les excitations indivi- 
duelles, il existe une énergie d’excitation e,. Seules 
les excitations du type spinon sont associées à des 
énergies faibles. Toutefois, ces spinons ne semblent 
pas susceptibles de transporter du courant ; 
d’ailleurs, statistiquement, le nombre d’excitations 
du type spinon est infiniment petit par rapport aux 
autres types d’excitations. Notre système est donc 
isolant. 

Ajoutons une dernière remarque au sujet des 
Spinons ; comme nous l'avons vu, les termes 
d'interaction proprement magnétique entre les ions 
voisins ont été négligés. Dans un véritable antifer- 
romagnétique au zéro absolu, l'énergie réelle des 
spinons serait donc assez différente de celle que 
nous calculerons dans la suite de notre exposé ; 
cependant la structure des excitations resterait de 
toute façon inchangée. 

Notre but étant défini, les pages qui vont suivre, 
seront consacrées aux développements mathé- 
matiques qui justifient nos assertions. 


d’excitations individuelles (5 — 1) et une exci- 


2. Excitations longitudinales et transversales de 


l’hamiltonien diagonal. Relations de commutation. | 


__ Tant que l’on ne tient compte que de la partie 
diagonale H, de l’hamiltonien, les excitations élé- 


mentaires gardent une forme simple, et le tableau | 


des opérateurs 5,4 et Su associés à la création ou à 
l’annihilation de ces excitations, a été dressé au 
paragraphe ITI-7. 


La question se complique quand on étudie le sys- | 


tème régi par l'hamiltonien À = H3, + Hypy. En 


effet, l'opérateur Hypy introduit des couplages | 


entre les différents opérateurs S et S. En pre- 
mière approximation, il est légitime de considérer 
les opérateurs de création et d’annihilation des 
excitations de l’hamiltonien H comme des combi- 
naisons linéaires et homogènes des opérateurs S% 
et Sue 

Avant tout, nous devons un peu préciser les pro- 
priétés de ces opérateurs. Ils appartiennent à deux 
groupes : 


a) OPÉRATEURS LONGITUDINAUX. — Îls sont de 
la forme Aÿ° B£ ou Bf° A. Nous poserons : 


Pp DE Lente Bi AE + eùk AE, BŸ] (57) 


lee - ae 
lx,k+0 ini AR B$+a + e—ik Bi Af+al. 


Ainsi pouvant être définis tous les opérateurs 
longitudinaux et de manière unique, à condition 
toutefois que k soit un vecteur de la zone réduite, 
tandis que le vecteur q peut varier dans toute 
l’étendue de la zon e de Brillouin. 

Les relations de commutation de ces opérateurs 
sont alors faciles à calculer, à condition de les 
évaluer en valeur moyenne pour un système à la 
température T.. 

Avec les conditions 4 € (DZ;) et q € (ZB), nous 
obtiendrons les relations approchées : 


Ms ia = 0 [sr + a lu,u + s] = 0 


Lire url = Op Vaslfe — 8k+al- (58) 


Les probabilités d'occupation des niveaux f#(T) 
et gk(T) ont été définies au chapitre IV-1. 


b) OPÉRATEURS TRANSVERSAUX. — Ils sont de 


la forme 4° A£ ou B+° Be (e et e’ représentant 
des états de spins opposés). 
Nous poserons : 


de AC PI ’ NE > 

trkta = Ar+a A% tar — en BE BE (59) 
+o@ 22e a À 7 

kk+a — Ax Ab+ r+g,k — ei Br Bêre 


Ces opérateurs sont définis de manière unique à 
condition que les vecteurs k et qg appartiennent à la 
zone de Brillouin réduite (DZ). 
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Les commutateurs de ces opérateurs sont très 
faciles à calculer et nous pourrons encore les expri- 
mer en valeur moyenne à la température 7. 

_ Avec les conditions ke(DZ;) et qe(DZ;), nous 
obtiendrons les relations approchées : 
+ + 
tte one = 0 
+ 
ler tuuts) = dku daslfk — 8k +] 
lo + ‘ + 4 La + 
Re ture] = 0 NA tuuts] = 0 (60) 
LA L + 
Cn+c, lunel = dyu Vostfr — Sk+ a) 

Par ailleurs, considérons les trois groupes d’opé- 
rateurs t et t*,t' et tt, L'et L*. Tous les commu- 
tateurs formés à partir de deux opérateurs appar- 
tenant à deux groupes différents sont nuls en 
valeur moyenne. 

Toutes ces relations nous servirent dans la suite, 
et déjà nous remarquons que les relations de com- 


(RG autel 0 


mutation n’introduisent pas de couplage entre opé- 


rateurs longitudinaux et transversaux. De même, 


nous verrons que dans l’hamiltonien, les termes 


longitudinaux et les termes transversaux peuvent, 
en un certain sens, être séparés. En conséquence, 
dans le système régi par l’hamiltonien 

H = Hp + Hypv, 


il existera encore des excitations longitudinales et 


._ . transversales représentées respectivement par les 


symboles (L) et (T). 

L'opérateur L associé à une excitation longi- 
tudinale et fonction linéaire des opérateurs l* et / 
sera choisi de telle façon que : 

[H, I] = Er L. (61) 

De même, l'opérateur T associé à une excitation 


. transversale et fonction linéaire des opérateurs 1* 


et t’ ou t'* et { sera choisi de telle façon que 
(HN = Er T. (62) 
Dans ces égalités, £; et Er représentent les éner- 
gies correspondant à la création ou à l’annihilation 
d’une excitation. 


3. Caleul de l’hamiltonien non diagonal utile 
H ypy.— Examinons maintenant la partie non dia- 
gonale de l’hamiltonien H et voyons quels termes 


nous devons retenir. Le critère est le suivant : un 


terme est utile si son commutateur avec l’un des 
opérateurs Si ou 5,4 redonne approximativement 
un terme analogue. 

Il est facile de voir alors que, de l’hamiltonien de 
Bloch H3y» (voir formule (7)), on ne peut retenir 
aucun terme intéressant. Tous les termes utiles pro- 
viennent de l’hamiltonien d'interaction A, qu'on 
peut exprimer sous la forme suivante (voir for- 
mule (47)) : 


1 : 
HE 3x2 g(g) (63) 


. € : P+E 
[cos Déto AE — sin Dégge—Ak+) Brfol 
: ; 4n 
[cos De DRE — SIN DE, e—iA(u—q) By—a 


. . £ £ . Sr € 
[cos D? 42 sin Peu BA][cos Dé Af— sin Péeixk BE]. 
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Vu la forme des excitations longitudinales et 
transversales, il est clair que nous pouvons extraire 
de A; deux hamiltoniens pr et Hiwpr : Hrxpz est 
bilinéaire par rapport aux excitations longitu- 
dinales'et Hzypr bilinéaire par rapport aux exci- 
tations transversales. 

Tous les termes de l’hamiltonien utile Hypy 
appartiennent soit à l’hamiltonien H;ypr soit à 
lP’hamiltonien H;ypy, mais certains termes de H;pv 
figurent à la fois dans H;ypx et Hinpr, et par consé- 
quent : 


Hipv Æ Hinnz + Hinpr. 


Il n’y a pas de termes de couplage entre exci- 
tations longitudinales et excitations transversales ; 
par ailleurs, les opérateurs associés à ces excitations 
commutent. En conséquence, les excitations longi- 
tudinales et transversales seront tout à fait sépa- 
rées, ce que l’on pouvait d’ailleurs prévoir a priori, 
puisque la composante 6, du spin total commute 
avec l’hamiltonien transformé. 

Si donc nous nous intéressons aux excitations 
longitudinales, nous pouvons toujours poser : 


Hp + Hnoo HD Hip. (64) 


Si, au contraire, nous nous intéressons aux eXCI- 
tations transversales, nous aurons : 


H => Hp + Hnpv © Hp + Hinpr. (65) 


Nous allons maintenant préciser la forme de ces 
opérateurs. 

La forme de l’hamiltonien transversal {}ypr est 
relativement simple et présente un caractère 
d'échange. Nous remplacerons donc partout 
g(u — k) par sa valeur moyenne g, comme nous 
Pavons fait dans le calcul de l’hamiltonien diagonal. 
En changeant k en k + p, u en uw + p et q en 
q + p,nous pouvons réduire les termes semblables. 

On en tire alors immédiatement : 


g à ë ’ 
Hinpr = ere Sn Den + Obllere orel 


sin (Dy+a + Du) A hou) 


me 2’ cos (D (D ) Lui ue :] 
6 N Ê ke + a k k,k+6 k +a,.k 


COS (Dur q — Du) per +, el 


avec les conditions ke(DZ,), ue(DZ,) et qe(DZ;), 
les termes diagonaux étant exclus. 

L’hamiltonien longitudinal H;ynz Comprend 
deux termes, un terme coulombien H;yprce et un 
terme d'échange H;vprx (seul le terme coulombien 
est pris en considération par Gell-Mann et Brueck- 
ner dans leur théorie des électrons libres, parce 
qu’il est seul important pour de fortes densités). 

Posons qg' — q +p, et dans les hamiltoniens 


48 
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Hinpce et Hinpze véduisons les termes semblables. 
Écrivons alors ces expressions sous une forme 
compacte, en étendant le domaine de variation du 
vecteur g à toute la zone de Brillouin (le symbole Z” 
sera utilisé pour dénoter ce type de sommation). 


Hinouc = + 5° g(g) sin (Ou ra — De) L'Érta — tax] 
sin (Dee D) ler e ee Lol 

HinD1z = — . 2” sin (+0 — x) Lera — tax] 
sin (Du+a — Du) Laure — ral 


Z” cos (Dr+o + x) Lra + lktaxl 


mt 
2N 
cos (Dur + Du) DORE SE Fépirl (67) 


avec les conditions ke(DZ.), ue(DZ,) et qe(ZB), 
les termes diagonaux étant exclus. 

Finalement on obtient : 

Hinpz = Hinpze + Hinpze (68) 

4. Excitations transversales. Théorie plus com- 
plète. — Examinons maintenant la forme exacte 
des excitations transversales (avec renversement de 
spin). La composante ©, du spin total étant une 
constante du système, il s'ensuit qu'il existe deux 
familles d’excitations transversales conjuguées 
l’une de l’autre et qui correspondent à Ao, = + 1. 

Nous poserons : 


+ ‘ 
Ta = Ex tétte E UE ral 


L 1e (69) 
Ta — > Xe Uek+e À Uk ral. 


Les coefficients ?, et u, devront être choisis de 
telle manière que les opérateurs représentent bien 
des excitations du système, c’est-à-dire les équa- 
tions : 

ATV EE POLAR TEE, (70) 
devront être identiquement satisfaites. 

Les opérateurs T et T' représentent alterna- 
tivement des opérateurs de création ou d’anmihi- 
lation ; tout dépend du signe de Æ. Si un opérateur 
de création appartient au groupe 7, son conjugué 
est l’opérateur d’annihilation correspondant et 
appartient au groupe 7”. La réciproque est naturel- 
lement vraie. 

Comme nous l’avons vu, nous devrons prendre 
comme hamiltonien, en première approximation : 


H © Hp + Hinpr 


(voir formules (21) et (66)). 
Le calcul des commutateurs de H et des divers 
opérateurs {* et £, 1'* et {est immédiat si l’on se 
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reporte aux relations de définition et de commu- 
tation de ces opérateurs (voir formules (49) et (50)). 


t+ 
k,k+a 


(Hp, tuto] = (exo ax ex) 
‘ fa 
Ê (Hp, ttktal = (exte te) 6e 
‘ 
Un, G+axl = + (CP ES ex) Ue+a,k 


(Hp, ok] — (+0 == ex) k+ak (71) 


CHinpr, tè+al = —# (fe — gr+a) Sin (Dk+a + Lx) 
Z'sin (dure + Du léure E turaul 


—+ (fr oi &k+ a) cos (Dx+a — D) 
Z' cos (Dur — Du) EE “x tel 
CHinor, k+al = —£ (fe — gx +0) Sin (Drra + Dx) 


Z’ sin (Dura + Du) Féére SE Wutasl 


q 
Le] 


(fe — ge + 9) COS (Dx+a — Dx) 


Z' cos (Dy +4 — Du) Cha — tuto] 
(Hinpr, tax) = (fr+a — 8x) Sin (Dxta + Dx) 
Z' sin (Dyrip + Du) Lire a te aol 
— À (feu — 87) COS (D+4 — Dr) 
Z'c0s (Dia 0,) te Rral 
[Hinpr, tra = u. e+a — 8x) Sin (Drya + Dx) 


Z’sin (Dyra + Du) PAR + tutoul 
—# (Peta — 8x) COS (Dr +4 — Dx) 


E' cos (Dusa — Du) légitra — turaul. (72) 


Remarquons toutefois que les relations de com- 
mutations dans lesquelles figurent H;ynr n’ont de 
sens qu’en valeur moyenne et que les commu- 
tateurs sont des fonctions de la température, par 
l'intermédiaire des probabilités f,(T) et g(T). 

Remplaçons maintenant dans les relations (70) 
les opérateurs T, et T, par leurs expressions déve- 
loppées. Les équations couplées qui déterminent à 
la fois les valeurs des énergies d’excitation E, et les 
coefficients A9 et y (à un facteur de normalisation 
près) en découlent immédiatement, puisque les rela- 
tions (79) doivent être identiquement satisfaites 
(dans la suite, nous écrirons 4, uz et E, et non Àgar 
Ur et EE, pour simplifier les notations). : 

Nous poserons : 


14 
o 
G — 


=R Z’ sin (D; + D) 


Le — 840) A — (fra — 8x) tx] 
Tee r. L' cos (D — x) 


Le — 82 40) Xe + (fe+a — ge) wl. (73) 


NS FX 


en 
DS) 


LT vo 


 N°8-9) 


dut RES  d Den ASUS CARE COR Dr TEA re es nb FER: FT ST ne TE 
É - (LE Es 2) PAR IE D à APTE ER hs 


Les coefficients À, et 4 s'expriment en fonction 
de cet +. 


Re 6 sin (Dr + 9 + D) + + cos (Dx+a — D) 
: ki — SN D) 


o Sin (Drsg + Dx) — + cos (Dr — D) 
Ex+g — € — E 


(74) 


Le report de ces expressions dans les équa- 
tions (73) donne deux équations couplées en cet 7, 
linéaires et homogènes. 

Posons : 


Fil£) 


I 


RN 2’ cos? (D +9 — Lx) 


[ fr — gk+a 


fr + a rn | 
Ek+a — x —- E 


+a — ke — E 


GE) = 1—£ sin? (De+q + D) 


KT En UE Rio 8e 
Ek+g — ER —ZE  Eprg — % —E 
HalE) = É E' sin (Or+a + D) COS (De ra — D) 
Îk — Ska Îr+a — 8x | 
DRE RAT IANE FRE — |. (75 
ER eee (5) 


Les équations en cet ts’écrivent : 
6 GalE) — + HafE) = 0 

— 6 H(E) + + F(E) = 0. (76) 
Ce qui conduit à l’équation aux valeurs propres : 
O(E) = FilE) Gi) — HE) = 0. (77) 
Nous pouvons un peu simplifier l'expression 
de F$(E), GE) et H}(E). En effet, dans la moitié 
des termes, nous pouvons changer k + q en —k 
et ken — (k + q). Si l’on tient compte alors de la 


symétrie des fonctions fx, g et ®, par rapport au 
changement k — — k, on peut écrire : 


N 
eh+a — Ex 
2 {D ® = 7 À 
COS? (Dx+ a x) E 
FOREST Er 7) 
AO N 2 k 5k 14 


: Ek+a —— Ex 
sin? (Da + x) Le —(ek+a — Fil 


HQlE) = 2 É Ex — ge +a) Sin (De+a + D) 


E 
cos (Drya — Dx) fe — (cha — nl: qi 


Les fonctions étant écrites sous cette forme, il est 
alors visible que la fonction O(E)est symétrique par 
rapport à £. À la racine Æ de l’équation O(E) — 0 
correspond toujours la racine — Æ comme prévu. 
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Aux valeurs positives de E correspondent des opé- 
rateurs de création, aux valeurs négatives des opé- 
rateurs d’annihilation. 

Pour déterminer les solutions de l’équa- 
tion O,(E) —0, nous devons placer l’une par 
rapport à l’autre les courbes y, = FE), G(E) 
et ya —= H?(E). Les pôles doubles des deux fonc- 
tions ont même coefficient : O(E) ne contient donc 
que des pôles simples et change de signe à chaque 
pôle. Examinons alors le comportement de cette 
fonction au voisinage de ses pôles. Pour simplifier, 
nous supposerons que les valeurs de q sont très 
petites ; le cas où [g| = 0 sera étudié et nous extra- 
polerons les résultats. 

Jusqu'ici nous n’avons fait aucune hypothèse 
sur les valeurs de +4. Pour étudier plus simplement 
l'équation O,(£) — 0, nous pourrons supposer que 
les fonctions +, sont choisies de manière à mini- 
miser les termes d’énergie libre qui proviennent de 
l’hamiltonien diagonal H, (formule (33) : 


tg 2ox(T) = Sd (79) 
SLR LE eo(T) 
PEUT NET 2 ox 


En conséquence : 


de FA 6 eo Sin 2@%x | 
LE) 1 22 . Z'(fx — 8x) ee AE PR 
D €o Sin? 2x | 
HE) = 1 +22 (fe 2 80 
C o(Æ) I 2 N (fu gx) ( sin 2 2 © LE eè ( ) 


Sin E sin? 2% 
BE) = FE DRl 
Remarquons, en passant, que ©,(0) = 0, car 
F${0) = H$(0) — 0 (cf. formule (32)). Cette racine 
correspond au spinon de vecteur d’onde nul : son 
énergie est nulle comme prévu. Soit alors : 


g eo Sin 2% 
MES 2 a , OT: 2 - A 
: N VE 8x) E? sin? 2 pr — ci 
, ; €o 
Lorsque Æ tend vers l’une des valeurs +————, 
sin 29% 


4 tend vers l'infini ; par ailleurs, au voisinage de 
ces valeurs, les fonctions F$(E), GQ(E) et H,{E) 
peuvent être remplacées par leurs équivalents. 


PFolE) © Fr + 
GolE) © Gr + sin? 2 pr tx 
H,(E) © Hg + sin 2 tr. 
On en tire alors aisément : 
O,(E) © (Fr + x) (Gr + sin? 2 9x tx) — (Hx + sin 2 9x tx)? 
te (sin? 2x Fr + Gr — 2 sin 2x Hx) = 


 T 


sin 2 @x 
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Ainsi se.trouve précisée la nature des pôles de la 
fonction O,(Æ). 

Par ailleurs, la forme de l’équation O,(E) — 0 
montre qu’il existe autant de racines que O,(Æ) 
possède des poles, et ces racines doivent être 
toutes réelles, puisque l’introduction d’un terme de 
couplage supplémentaire Hypy doit laisser le 
nombre d’excitations inchangé. Cette remarque 
permet alors de tracer la courbe y = O,(E) (voir 


fig. 10). 
1@, (E) | 


1 


=, 


Excitation 
individuelle 


| 


Fic. 10. — Excitations transversales. 
Représentation de OQ(E). 


Nous observons alors que toutes les racines de 
l'équation O,(Æ) — 0 sont voisines des pôles 
de O,(E), à l'exception des deux racines, les plus 
basses en valeurs absolues qui, pour de faibles 
valeurs de |g, sont nettement séparées de l’en- 
semble continu des autres racines ; ces racines is0- 
lées correspondent à une excitation de caractère 
collectif, comme on peut s’en rendre compte faci- 
lement, si l’on calcule les coefficients ?4 et uy qui 
lui sont associés : ce sont des spinons. 

Remarquons alors que, pour que ces spinons 
existent réellement, il faut obligatoirement que 
O,(0) > 0, au moins pour de petites valeurs de |q{|. 
Cette condition est satisfaite si les fonctions o4 
obéissent aux égalités (79), mais il n’en serait pas 
de même si l’on donnait aux fonctions œ des 
valeurs arbitraires. 

En résumé, les énergies des excitations indi- 
viduelles qui forment un ensemble quasi continu, 
ne sont que faiblement déplacées, et la limite infé- 
rieure de ces énergies d’excitation, pour de petites 
valeurs de g, est précisément e,. Il existe, par 
ailleurs, une excitation collective ou spinon associée 
à chaque vecteur q. 

L'énergie du spinon croît à partir de zéro 
comme |q|, puis, au delà de certaines valeurs de q, 
se fond dans l’ensemble continu des excitations 
individuelles. 

L'énergie du spinon, pour de faibles valeurs 
de |q|, peut être facilement calculée, et, pour T —0, 
l’on obtient : 

E — qC(q19) 
avec 


5 Egg) PUITS 
ip | VUE PRES, ee 2 
C’{n) Fe one pe 2’ Sin 2quir VU) | 


2g ; 2g ; ei 
AS sin Lor cos Doll ss 90 
Res sin 20 Cos 2 qe. he sin 2 qe. 61) 
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5, Existence des spinons et validité de la théorie. 
__ Comme nous l’avons vu, pour que toutes les | 
excitations aient un sens, il faut que @,(0) > 0; | 
et en particulier que @,(0) >, 0 ; nous avons 
vérifié par ailleurs que cette dernière condition 
était satisfaite dans le cas où les +4 obéissaient 
aux relations (79) qui correspondent à la mini- 
misation de l'énergie libre dans le cas où l’on néglige 
tous les termes non diagonaux de l’hamiltonien. 

Abandonnons cette restriction et donnons aux 4 
des valeurs arbitraires mais cependant telles que la 
condition @,(0) > 0 soit remplie. Nous pouvons 
alors envisager de calculer les corrections à apporter 
à l'énergie libre, du fait de l'introduction des 
termes Hypy ; la méthode utilisée serait analogue à 
celle qu'ont développé Brueckner, Sawada et Gell- 
Mann pour un gaz d'électrons libres. L'énergie 
libre corrigée s'exprime alors en fonction des oz qui 
peuvent être alors déterminés en minimisant cette 
énergie libre. 

Or la condition @,(0) > 0 restreint les variations 
de ®,. Par conséquence, l’énergie libre du système 
possède des extrêma sur la surface O,(0) — Otracée : 
dans l’espace des ®,. Il semble alors légitime 
d'admettre que le minimum de l'énergie libre 
correspond effectivement à l’un des points de cette 
surface. Cette condition traduit le fait que, lorsque 
le vecteur |g| tend vers zéro, l'énergie des spinons 
tend vers zéro, ce qui semble raisonnable. 

Examinons alors les conséquences de cette 
remarque. Soit ow les valeurs de vw, qui satisfont 
aux conditions tg 20% — €pler (e9 étant lui-même 
fonctions des wz). 

Admettons que l'introduction des termes Hypv 
crée seulement une légère perturbation, les valeurs 
des +, seront voisines des valeurs ow. Or, pour 
T — 0, on peut écrire en toute généralité : 


@5(0) — F6(0) Go(0) 
et de plus 


[Go(0)le,=e,x N 


8 5 sin? 2.001 > "0: 
€o 


Donc, si l’on fait varier les valeurs de oz au 
voisinage des valeurs o4, la fonction G,(0) reste 
positive. La condition @,(0) — 0 s’écrit donc: 

28 ” 4 


RO) == 
à Es N à = ÿ 


avec 
% — Ex — 6x COS 2x + epsin 2x (formules (22) et (27) 
et 


Nous en tirons la relation : 


Es] DRASS ES des 
| dpx 9x Por Œ Fe ei 2qu] Se 
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et, par conséquent : 


(ous. a Z’ sin? 2 ou des. (82) 
. Si, par conséquent, l'introduction de termes non 
diagonaux peut faire varier légèrement les para- 
mètres 7, l'énergie d'activation e, doit toutefois 
rester inchangée, puisque nous avons des raisons 
de croire que la condition @,(0) — 0 doit rester 


M tistaite. 


RU ct il n’en est pas de même lorsque 
le coefficient g devient petit par rapport à la 
largeur de la bande de Bloch. Dans ce cas, si l’on 
tient compte des termes non diagonaux, l énergie 
associée à un état de base métallique peut devenir 
plus basse que l’énergie de l’état de base antiferro- 
magnétique ; il y à alors, dans ce cas, transition 
brusque de l’état isolant à l’état métallique. 


6. Excitations longitudinales. Théorie plus com- 
plète. — Examinons maintenant la forme exacte 
des excitations longitudinales (avec conservation 


- duspin). Nous poserons : 


Lis lholente ue bia (83) 


Le vecteur q est un vecteur quelconque de la 
zone de Brillouin. Les coefficients à4 et u; seront 
choisis de telle manière que les opérateurs Z, repré- 
sentent bien des excitations du système, c’est-à- 
dire que les équations 


Ed ER (84) 


devront être identiquement satisfaites. 

Les opérateurs L représentent alternativement 
des opérateurs de création ou d’annihilation, tout 
dépend du signe de Æ. Les opérateurs de création 
ou d’annihilation sont naturellement conjugués 
comme on le vérifie sur l’équation (84). 

Comme hamiltonien, nous devons prendre, en 
prernière approximation : 


HIDE INDE 1) et (67)). 


Le calcul des commutateurs de Æ et des opéra- 
teurs L* et L est immédiat, si l’on se reporte aux 
relations de définition et de commutation de ces 
opérateurs (voir formules (57) et (58)). 


(voir formules (2 


Cp, Rrta) = (ex+a he) let 


(Hp, Riakr = (Ek+a — €&%) +oa 


2 e(q) — 
CHinot, léxta] — pets] (x — gr + a) Sin (Dr +4 — Dr) 


Z'sin (Dyra — d ere -z ly+a, “] 


_. £ (fe — gk+a) COS (Dkra + Dx) 


Z' cos (Quyo + Du) Lure + lutuu] 


TRANSITIONS ENTRE ÉTATS MÉTALLIQUES ET ISOLANTS RTE 


tt Ü+ a — 8) Sin (D + à — D) 


lutau] 
+ À (he — 84 a) COS (Déra + De) 


(HinDz, ra) = | 


no, Ti 


2’ cos (Dyra + Du) Lhuta + lutqu]. (85) 


Remarquons que les relations de commutation 
dans lesquelles figurent H;yp; n’ont de sens qu’en 
valeur moyenne et que les commutateurs sont des 
fonctions de la ni par l’intermédiaire des 
probabilités f,(T) et g,(T). 

Remplaçons ae dans la relation (84) 
l’opérateur L par son expression développée. Nous 
en tirons des équations couplées qui permettent à la 
fois de calculer les valeurs de l’énergie d’inter- 
action Æ£ et des coefficients À; et uy. 

Nous poserons : 


2g(q) — g1 &, . 
5 = [ETS] E'sin (Qui o — Eu) 
[fe — Su+ a) lu + 


T = F2 COS (Dur a + Du) Lu —gu +a) Me — 


(fu+ a — Su) Lu] 


(fu + à — 8u)u] 
(86) 


Les coefficients 24 et u4 s'expriment en fonction 
de cet de +. 


— 6 Sin (Dr — D) + + cos (Der + D) 
cktq x — E 
6 Sin (Dr+a — Dr) + + cos (Drro + D). 
CET Dome €, —E 


À = (87) 


Uk = 


Le report de ces expressions dans les équa- 
tions (86) donne deux équations couplées en cet 7, 
Les équations en cet rs’écrivent : 


5 FE) —< (EE) mur) 0 (88 


o HAE) £r GE) = 0. 
Ce qui conduit à l'équation aux valeurs propres : 


GE) = FE) GE) + [OA] Hàte) = 0. (89) 


Après simplification les expressions F{(E), G(E) 
et H,(E) prennent les formes suivantes : 


RUE) = 12 (EE) je — sd) 


sin? (Dx+g — D | Es ï | 
Fra i E? — (cx+g — x)” 
GalE) = 1 + 2 É D'fe — gr+0) 


= LA 
Ek+a — € 
cos? (Dig + )| ts. | 
( L+a fc SRE por x)? 


Z’ sin (Dxra — Lx) 


Ek+g —— € 
e Pr 10 [ss | (90) 
cos (Dr+a + D) Fine le 
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Les fonctions étant écrites sous cette forme, il 
est visible que la fonction O(Æ) est symétrique par 
rapport à Æ. Aux valeurs positives de Æ corres- 
pondent des opérateurs de création d’excitation, 
aux valeurs opposées de Æ des opérateurs d’anni- 
hilation qui leur sont associés. 

Remarquons encore que @,(E) ne contient que 
des pôles simples, et par conséquent, change de 
signe à chaque pôle. Le nombre de racines de l’équa- 
tion O,(Æ) — 0 sera encore égal au nombre de pôles 
de la fonction O,(E). 

L'étude générale des excitations étant difficile, 
étudions simplement le cas où q — 0. Nous suppo- 
serons que les fonctions w, obéissent aux éga- 
lités (78) et que la température Test nulle. 

Dans ce cas simple, nous pourrons poser : 


METIER 
Ci g 
Les différentes fonctions s'expriment alors sim- 
plement en fonction des +4. 


DM D NUE (91) 


n 


= !/ er DA ®P.. 
1 N (vY D) 


à eo Sin 2 


F,(E ——— : 
o(E) 13 sin? 29x 2 € 


& +, Co Sin 2x COS? 29% 


HE) = 1 +2É2 4 GE 
C o(Z ) 1 } 2 N EL? sin? 20% be eÿ | 
1 4 g eo Sin 26p COS 20% 

E G RE DEN ®D.. S 0. 
xt], NOV 0 Er sine 29 — eà 


L'équation aux valeurs propres s'écrit alors : 
Oo(E) = Fo(E) GolE) = 0. 


Les valeurs propres Æ£ sont donc solutions de 
June ou l’autre des équations F$(E) — 0 ou 
GE) = 0. Dans chaque équation, le nombre de 
racines est égal au nombre de pôles des fonctions 
F(E) et GE), mais les fonctions F,(E) et G,(E) 
comportent moitié moins de pôles que les fonctions 
FE) et G(E), puisque, pour qg = 0, nous aurons 


LA La 
Et+q —— EE = E—k+eg — E—X- 


Le nombre total de racines est donc bien conservé 
comme il se doit. 

Les allures des fonctions F$(E) et G,(E) sont 
représentées sur les figures 11 et 12. 

Examinons les racines de l’équation G,(E) — 0. 
Nous remarquons que Gi(es) = 0, et que toutes 
les excitations qui correspondent à ces racines ont 
le caractère d’excitations individuelles, caracté- 
risées par une énergie d'activation e, et forment un 
ensemble quasi-continu. 

Étudions alors les racines de l'équation 
F(E) = 0. À ces racines correspondent des exci- 
tations individuelles qui forment un ensemble 
quasi-continu et une excitation collective. Comme 
toujours, pour les excitations individuelles, il 
existe une énergie d'activation e,. l'excitation 
collective, quand elle existe, est du type plasmon, 
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Son énergie est nettement supérieure à celle de 
toutes les excitations individuelles. 


1 


Fic. 11. — Excitations longitudinales. 
Représentations de FG(E). 


F1G. 12. — Excitations longitudinales. 
Représentation de GE). 


L'énergie du plasmon peut se déduire de l’équa- 
tion implicite en Æ (formules (91 et 92)) : 


2e? Je eo Sin 20% ; 
®D,)? L GES IS 93 
TND Z (VV D) E? sin? 2x — e$ : (93) 
Nous pouvons évaluer approximativement 


l'énergie £ du plasmon en procédant comme au 
paragraphe IIT-6 et en définissant une énergie de 
coupure €. 

La formule (105) se transforme et donne : 


PR ee pe AA 
0 


GT n E2— eÿ —e? 


de = 1 (94) 


(S représente la surface de la zone réduite dans 
l’espace des vecteurs d'onde). 

En particulier, dans le cas où E£? > e? + €, on 
obtient approximativement : 


2 2 2 
\Veg + eè 


L'énergie des plasmons décroit donc avec @& 
comme on pouvait s’y attendre. 

En résumé, nous observons que, pour toutes les 
excitations individuelles, l'introduction dans l’ha- 
miltonien des termes non diagonaux Hmpy laisse 
subsister une énergie d’activation eo. Toutefois, à 
côté des excitations individuelles, il existe des 
excitations de caractère collectif, spinons et plas- 
mons, Inais ces dernières ne détruisent pas le caräc- 
tère isolant du système. 


* 


\ 


a 


RCE 


De | 
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VIII 


APPLICATIONS DE LA THÉORIE 
AUX MODÈLES DÉSORDONNÉS 


1. Remarques générales. — Jusqu'ici nous avons 
supposé que le réseau ionique qui constitue notre 
modèle est ordonné et facilement décomposable en 
deux sous-réseaux également ordonnés. Nous avons, 
de plus, admis que le nombre des électrons était 
égal au nombre des ions. 

Or, dans les antiferromagnétiques non stoechio- 
métriques comme Ti20* (n # 3), les atomes cons- 
tituants ne peuvent former un réseau ordonné. 
Dans les bandes d’impureté, il en est de même, et 
les centres d’impureté sont en général répartis tout 
à fait au hasard. De plus, les spécimens semi-con- 
ducteurs contenant le plus souvent plusieurs types 
d’impuretés, il y a compensation partielle. Ainsi, 
par exemple, aux basses températures, dans un 
spécimen de type n, plusieurs électrons provenant 
des niveaux donneurs sont captés par les impuretés 
acceptrices. Le nombre des électrons qui se dépla- 
cent dans la bande formée par les impuretés dona- 
trices est donc inférieur au nombre des ions qui 
constituent le support de cette bande. L'expérience 
montre alors que, même dans ce cas, tout se passe 
comme s’il existait une énergie d’excitation pour 
les électrons de la bande d’impureté. 

Ce phénomène est difficile à expliquer si l’on part 
d’un schéma de bandes rigides. Un essai d’inter- 
prétation a toutefois ‘été donné par Mott [5]. Il 
part du fait qu’un centre minoritaire chargé doit 
retenir autour de lui une charge égale et opposée de 
façon à sauvegarder la neutralité électrique. Ainsi 


- pour un spécimen de type n, si un centre accepteur 


capte un électron, il devient chargé, et doit retenir 
un trou. Il admet alors que ce trou ne peut devenir 
tout à fait libre que si on lui fournit une énergie 
suffisante pour l’arracher au centre accepteur. 
Cette interprétation, toutefois, semble douteuse : 
on peut, en effet, admettre que la charge compen- 
satrice est simplement due à la polarisation des 
électrons de la bande d’impureté due aux centres 
donneurs. 

Reprenant alors la question sous un angle nou- 
veau, nous verrons que, par application directe de 
notre théorie au cas désordonné, on peut montrer 
que, pour de faibles densités en impureté, les élec- 
trons de la bande d’impureté ne peuvent être 
excités sans énergie d'activation ; c’est cette ques- 
tion que nous allons maintenant étudier. 


2. Extension de la théorie. — Soit alors un 
grand nombre d’électrons en mouvement dans un 


, 
à 


potentiel aléatoire, créé par un réseau désordonné 
d'ions et par les charges moyennes dues aux élec- 
trons, Les solutions du problème à un électron sont 
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des ondes analogues aux ondes de Bloch. On peut 
encore définir approximativement pour chacune de 
ces ondes un vecteur d'onde k. Dans le cas à une 
dimension, la correspondance entre un vecteur k 
et sa fonction d’onde est automatique, même si le 
potentiel n’est pas périodique ; en effet, le nombre # 
peut être relié au nombre de zéros des fonctions 
d’onde supposées réelles. Dans Je cas tridimen- 
sionnel, la situation est plus délicate ; en fait, la 
plupart des méthodes d’approximation utilisées 
pour calculer les fonctions d'onde dans l’espace 
à trois dimensions, font appel plus ou moins impli- 
citement à la notion de vecteur d'onde. Dans le cas 
général, on peut donc, sans grande erreur, associer 
à chaque onde son vecteur d'onde k. 

Simplement, l'introduction d’un potentiel aléa- 
toire tend à élargir la bande de Bloch en favorisant 
la création de queues de bande vers les hautes et 
aussi les basses énergies. 

Formellement, l’hamiltonien de Bloch est donc 
identique que l’on se place dans le cas ordonné ou 
le cas désordonné. Toutefois, la notion de zone de 
Brillouin perd toute signification nette, et les 
niveaux E, ne sont plus périodiques par rapport à ke. 

L’hamiltonien d'interaction garde lui aussi une 
forme analogue, mais la conservation des vecteurs . 
d'onde n’est plus qu’approximative. Plus préci- 
sément, dans un réseau ordonné, l’hamiltonien 
d'interaction contient deux sortes de termes : les 
termes normaux pour lesquels la somme des vec- 
teurs d’onde appartiennent à la première zone de 
Brillouin est strictement conservée, et les termes 
«€ Umklapp » pour lesquels cette somme n’est con- 
servée qu'à une translation du réseau réciproque 
près. Si donc le réseau est perturbé, les termes nor- 
maux seront à peu près conservés, mais il est nor- 
mal de penser que les termes « Umklapp » seront 
très sensiblement modifiés. 

Voyons alors comment notre méthode peut être 
appliquée au cas désordonné. Soit S la surface de 
Fermi associée à l’énergie du type de Bloch et V le 
volume de l’espace réciproque intérieure à cette sur- 
face. La surface de Fermi, en général, n’est pas 
sphérique. En effet, un réseau plus ou moins désor- 
donné conserve souvent des directions cristallines à 
peu près définies comme cela se rencontre dans les 
antiferromagnétiques. De même, dans les semi- 
conducteurs où les centres d’impureté sont répartis 
complètement au hasard, la surface de Fermi asso- 
ciée à ces centres possède la symétrie du cristal 
porteur. La forme de la surface de Fermi pourra 
rappeler celle d’un polyèdre simple ; par exemple, 
dans le cas d’un cristal cubique, le polyèdre sera un 
cube ou un octaèdre. 

Ces polyèdres seront tels que les faces opposées 
puissent se correspondre deux à deux par trans- 
lation. Soit p, ces translations. Pour un cube, par 
exemple, nous aurons trois translations p;, Poe et Ps 
parallèles aux arêtes du cube. Le volume du cube 
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sera pris égal à V, c’est-à-dire qu’il sera tel qu’il 
contienne autant de niveaux qu’il y a d’électrons 
dans Ja bande que nous considérons. Le polyèdre" 
choisi jouera alors le rôle de zone réduite. 

* Dans le cas ordonné, à chaque vecteur k de la 
demi-zone réduite, était associé un vecteur & + p 
extérieur à cette zone. De même, ici, nous asso- 
cierons à chacun des vecteurs Æ intérieur au 
polyèdre, l’un des vecteurs k + p,. Contrairement 
au cas ordonné, les différents vecteurs p; ne sont 
pas équivalents entre eux, de même que les vec- 
teurs 2p, ne sont pas équivalents à zéro. Par consé- 
quent, le signe et l’indice 7 du vecteur choisi dépen- 
. dront de la région du polyèdre à laquelle appartient 

le vecteur k. Ainsi un cube de centre O sera décom- 

posé en 6 pyramides de centre 0 ayant pour bases 
les faces du cube. Les faces opposées du polyèdre se 
correspondent point par point. Le vecteur corres- 

pondant à k sera noté k’. 

Nous pouvons alors effectuer une transformation 
analogue à la transformation définie en IIT-3 (à 
condition de remplacer k + p par k'). Comme 
précédemment, nous obtenons ainsi un abaissement 
de l'énergie de base du système. Ce résultat est 
particulièrement facile à prévoir si l’on se réfère au 

concept de paires, comme au paragraphe VI. Le 
gain d’énergie provient de l’échange entre paires. 

Les termes directs et les termes € Umklapp » parti- 

cipent à cet échange. 

Dans le cas ordonné. le choix du vecteur p est 
impératif, parce qu’alors les termes directs et les 
termes « Umklapp » doivent apporter leur contri- 
bution, pour que l’énergie de l’état de base soit 
minimum, Il n’en est pas de même dans le cas désor- 
donné ; en effet, la conservation rigoureuse des 
moments n’est pas imposée aux termes «Umklapp », 
et la contribution de ces termes est très faible. Au 
contraire, le gain d'énergie obtenu grâce à la trans- 
formation est d'autant plus important que l’échange 
porte sur un plus grand nombre de paires. Il faut 
que la demi-bande inférieure soit pleine. C’est 
pourquoi le volume du polyèdre n’est pas arbi- 
traire ; il doit contenir autant de niveaux qu'il y a 
électrons dans la bande, Les valeurs des p; doivent 
donc s’ajuster automatiquement de façon que la 
demi-bande inférieure soit pleine, ce qui conduit à 
une minimisation de l’énergie de l’état de base. 

Ainsi l’on peut comprendre pourquoi dans les 
bandes d’impuretés aux faibles densités, une éner- 
gie d'activation d’impureté se manifeste toujours, 
même dans les spécimens partiellement compensés. 

Du point de vue pratique, l’état de base peut être 
calculé en ne tenant compte que des termes dia- 
gonaux de l’hamiltonien. La forme de AH, reste 
pratiquement inchangée. Toutefois, les termes 
d'échange sont scindés en plusieurs groupes, chaque 
groupe étant associé à l’un des vecteurs p,. Les 
vecteurs d'onde appartenant à la région j seront 
notés k; et leurs vecteurs conjugués seront 


JOURNAL DE PHYSIQUE | . 


N° 8-9 


k! —k, + p;. Les ‘termes € Umklapp » étant 


.négligés, les seuls échanges possibles seront du type 


normal (voir figure 13 et n’interviendront qu’entre 
paires appartenant à la même région 7). 


Fig. 13. — Échanges normaux pour un réseau désordonné. 
L'hamiltonien diagonal s'écrit : 
[æ) 
— Z'[(Ex cos? où + Ex sin? ox) Aë A 


+ (By sin? of + Er cos? et) Br ° Bt] + Hion 


S s, ES LE +e ne +e 4e Serre 
— 7 N > [Ax A% + Pr Bx] Bis Ay + Puy Bi] 
N y : 
œ 
DRASS NO ALES NE k) 
a D 2x (Ar At; — Br; PE, 


k.,u. 
ji 
sin 2 Qu (Au; Au; — Ba Be) | (96) 


(la notation Y indique que la sommation se fait 
kj 


J 
sur tous les vecteurs & intérieurs à la zone 7j). Le 
terme Hop est un terme coulombien sans impor- 
tance pratique dans notre calcul. 
La fonction d'onde de l’état de 
toujours : 


base sera 


PIN AS CES Me 


Tous les zones 7 étant supposées identiques, on 
pourra définir une énergie d’excitation e, unique : 


8 


È (97) 


© Ep = €; = a sin 2 ke 


J 


L'énergie de l’état de base sera minimum, si l’on 
a encore : 
avec Ex — Er ee Ey. 


€ > 
t8 2x =— (98) 


ke 

On remarque que l'énergie d’excitation est plus 
faible que dans le cas ordonné. Supposons, par 
exemple, le polyèdre cubique : tout se passe comme 
si le coefficient g avait été remplacé par g/6. 
L’omission des termes « Umklapp » en est la cause 
directe. 

Notre théorie s'applique done au cas désordonné 
moyennant quelques modifications ; ce résultat n’a 
rien d'étonnant. En effet, nous avons vu précé- 
demment que le gain d’énergie obtenu provient de 
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… la ségrégation partielle d'électrons de spins opposés 
- sur deux sous-réseaux différents. Or, pour que cette 
“ condition soit réalisée, peu importe que les réseaux 


[M, soient ordonnés ou non et peu importe que le 


… nombre des électrons soit vraiment égal au nombre 


… des ions, pourvu que la différence soit petite. 


É 8 
Toutefois, et comme dans le cas ordonné, pour 


— que ce processus soit effectivement réalisé, il faut 


… que l'énergie du système dans cet état de ségré- 
—… gation partiel soit inférieure à l’énergie de l’état 


| métallique, compte tenu de l'énergie de corrélation 
|, classique. Cette énergie de corrélation est celle que 


- l’on peut calculer, quand on tient compte de la 
- modification de l’état de base due à l'introduction 


| - des termes Hypy définis au paragraphe VIII. 


IX 


DENSITÉ DES IONS 
M ET PASSAGE DE L’ÉTAT MÉTALLIQUE 
A L’ÉTAT ANTIFERROMAGNÉTIQUE PARTIEL 


kr L'énergie de l’état de base du système supposé 
. désordonné ou non possède en général pour des 
» variations des fonctions +, plusieurs extrema. 
L'état métallique (ox == 0) correspond toujours à 


— l’une de ces valeurs extrêmes de l’énergie. Cette 


… propriété est évidente, puisque tout changement 


… des w, en — vw; laisse l'énergie de base invariante ; 


…_ en effet, les spins associés aux sous-réseaux À et B 

sont implicitement interchangés. Pour de fortes 
- densités ioniques, cet état métallique est stable ; 
c’est l’état de base du système. 

Par ailleurs, quand la densité décroît, il apparaît 
un état antiferromagnétique partiel tel que l'énergie 
de base associée à cet état soit aussi minimum pour 
de petites valeurs de vw. Une première condition 


” 
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électronique critique, il y a passage brusque de 
l’état métallique à l’état isolant. 

Pour nous fixer les idées, nous pouvons examiner 
le comportement d’un gaz d'électrons libres ; la 
densité est définie par le nombre sans dimension r, 
(voir I-1). Essayons alors d’estimer la valeur de r, 
qui correspond au point de transition. (Nous note- 
rons : do — h?/me? rayon de Bohr ; E, —#?/2ma2 
énergie d’ionisation de l’atome d'hydrogène). 

L’énergie des électrons libres a été calculée par 
Gell-Mann, Brueckner et Sawada [4]; on obtient 
pour l’énergie de base, par électron : 


Er, = Eg + Er + Ex + Ec (99) 
E£ est l’énergie de Bloch 
Eg 3/9r\2/5 4 2,21 
nn le . 
Ex est l’énergie d'échange 
Er 3 [9r\ 2/31 0,916 
nie Le 


Ex’ est aussi un terme d’échange, obtenu par un 
calcul de perturbation au deuxième ordre ; il vaut 


Er les ="0 046 (102) 
E est l'énergie de corrélation 
EG æ 3 CSS JEE du 
E; NA 1 GE —Co ètre 
a rs a rs 
| Log [1 + LE Quui] LE Qutui 
avec 
: +09 F : q° 
Qu= [dr f 7 aemmexs [1 (LT +ap)|. 
Ip] <1 
Ip+ql>1 


Les variations de la fonction Æ£.(r,) obtenues par 
intégration numérique sont représentées sur la 
figure 14. Remarquons que les approximations 


; pour qu’un tel état puisse apparaître, a été donnée 
…—_ au paragraphe IV-7 ; pour que notre modèle pré- 
… sente réellement une bande interdite de largeur n 
- (n < 6), il faut que l’énergie e, soit supérieure à la 
4 fluctuation maximum de l’énergie de Bloch sur la 
- surface de la zone réduite. Cette condition ne peut 
… être réalisée que si la densité électronique est assez 
4 basse, mais ce n’est pas un critère suffisant pour 
» affirmer que l’état antiferromagnétique obtenu pos- 
sède réellement une énergie inférieure à celle de 
l’état de base métallique. 

En effet, pour comparer les énergies des deux 


… états, il faut tenir compte de l’énergie de corré- 

“ Jation électronique. La réduction de l’énergie obte- 

: nue ainsi est due à un effet de polarisation des élec- 05 rs 
D trons et n’a été calculée que pour un gaz d'électrons 0 1, 2 8 4 5 6 7 8 3 10 

- libres. FiG. 14. — Énergies comparées d’un système d’électrons 


libres dans l’état métallique (£r) et dans l’état isolant 


Toutefois cette énergie de corrélation est certai- 
antiferromagnétique (Æ1). 


nement bien plus importante pour un état métal- 
“ Jlique que pour un état antiferromagnétique partiel. 
“ Il parait donc assez normal d'admettre avec 
- Mott [5] qu’au voisinage d’une certaine densité 


faites sont surtout valables pour de petites valeurs 
de r, ; néanmoins, par extrapolation, nous pouvons 
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utiliser la formule (99) pour estimer l’ordre de 
grandeur de la valeur de r, au point de transition. 
 Calculons maintenant l'énergie de l’état de base 
du système antiferromagnétique ; l’énergie de cor- 
rélation qui est faible, sera négligée. La première 
tâche consiste alors à évaluer le coefficient g. Nous 
remarquons sur la formule (21) que l’énergie 
d'échange vaut Er = — gJ4. 

Identifions cette valeur à la valeur de Æ, cal- 
culée pour un gaz d'électrons libres (101). On 
obtient ainsi une valeur très approximative de g, 
soit : 


15 (103) 


° Fe 
D 


6 Ee Es 3,66 


Tr \ 4 Te 7e 


Nous pouvons maintenant supposer que, pour 
un gaz d'électrons libres antiferromagnétiques, le 
réseau est cubique centré, ce qui implique que la 
zone réduite est un cube d’arête p. 

Posons : 


(104) 


L'énergie d’excitation e, et l’énergie de l’état de 
base Æ, sont données par les équations implicites 
suivantes : 


€ 2 Ê 
AU É (s: —;) Are sh (:) 


RES 
ie : PES 


1 15 3 o 
| £. (105) 


n 


Les variations de la fonction Æ,(r,) sont repré- 
sentées sur la figure 14. Le point de transition se 
détermine par l'équation ÆE,(r,) — Ej(r,). On 
trouve ainsi r, — 4. Cette valeur est évidemment 
trop petite, puisque une telle valeur correspond aux 
densités électroniques que l’on rencontre dans les 
alcalins ; toutefois, nous remarquons que la valeur 
de r, au point de transition dépend fortement de la 
valeur de g choisie, Pour déterminer le point de 
transition avec quelque précision, il faudrait cal- 
culer très soigneusement l’énergie d’échange. Il est 
toutefois vraisemblable que le point de transition 
corresponde aux valeurs de r, voisines de 7 ou 8. 


X 


ANTIFERROMAGNÉTIQUES 
ET BANDES D’IMPURETÉ. 
GÉNÉRALISATION DE LA THÉORIE. 


Les résultats obtenus sur les modèles simples 
étudiés Jusqu'ici permettent de mieux comprendre 
les anomalies de conductivité que l’on peut observer 
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dans les antiferromagnétiques et les semi-conduc- 
teurs impurs, phénomènes décrits dans les premiers 
chapitres. 

.Ces anomalies s'expliquent qualitativement, si 
l’on admet que, par suite de l'interaction coulom- 
bienne entre électrons au-dessous d’une certaine 
densité d’ions et d’électrons, il apparaît une faille 
dans les bandes d'énergie, la largeur. de la faille | 
dépendant à la fois des densités électroniques et de 
la température. 

Le mécanisme d’un tel phénomène a été analysé 
sur un modèle simplifié. Malheureusement, l’appli- 
cation directe de la théorie à des exemples concrets 
reste difficile. En effet, le problème se complique du 
fait que les états associés aux orbitales atomiques 
sont en général dégénérés ou presque dégénérés. 
Il faut donc admettre que les électrons peuvent 
occuper plusieurs orbitales sur le même ion et non 
plus une seule comme nous l’avions supposé 
jusqu’ici. | 

Cette situation se rencontre dans les antiferro- 
magnétiques. En effet, les niveaux intéressants 
appartiennent à la couche 3d ; il y a donc dégéné- 
rescence. En fait, comme nous l’avons vu, l’exis- 
tence d’un champ cristallin lève certaines dégéné- 
rescences ; néanmoins, au centre de la zone de 
Brillouin, il subsiste toujours des dégénérescences 
compatibles avec les symétries du cristal. 

Une situation analogue se présente dans les semi- 
conducteurs, comme le germanium et le silicium. 
Examinons d’abord le cas des semi-conducteurs de 
type n. La bande de conduction présentant tou- 
jours un certain nombre de vallées équivalentes, il 
s'ensuit que les niveaux d’impuretés donneurs 
doivent être fortement dégénérés. 

En fait, les fonctions d’ondes électroniques rela- 
tives aux centres d’impureté forment des représen- 
tations du groupe des symétries du cristal. 

On obtient donc ainsi, pour chaque niveau de 
type hydrogénoïde, autant de niveaux qu’il y a de 
vallées. Tous ces niveaux sont dégénérés ou presque 
dégénérés ; toutefois, il est vraisemblable que le 
niveau fondamental soit légèrement séparé des 
autres niveaux et corresponde aux valeurs &«;, — 1. 

Ainsi, seule la fonction la plus symétrique garde 
une valeur finie non nulle au voisinage du centre 
d'impureté. Or, à proximité de l’impureté, l’effet 
d'écran tend à diminuer ; l’électron qui gravite 
autour de l’impureté est done plus fortement attiré 
aux Courtes distances. L’amplitude de la fonction 
d’onde sur l’atome d’impureté est donc plus forte 
que ne le prévoit la théorie classique. Il en résulte 
donc un abaïissement de ce niveau, ce qui lève la 
dégénérescence. 

. À priori, cet effet doit être particulièrement sen- 
sible dans le silicium, car la dimension des orbites 
y est plus petite que dans le garmanium. De fait, 
pour le silicium, l’abaissement d'énergie est de 
l’ordre de 0,013 eV (voir Kohn [19]). 
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Cette énergie, toutefois, est encore de l’ordre de 
_ grandeur de l’énergie d’excitation des électrons qui 
se déplacent dans la bande d’impureté. En consé- 
quence, il est difficile de négliger les interactions 
entre les différents niveaux qui découlent de la 


= multiplicité des vallées de la bande de conduction, 


et particulièrement en ce qui concerne le germa- 
nium. La structure des bandes d’impureté, compte 
… tenu des interactions coulombiennes, doit donc être 
ù assez compliquée pour de faibles densités de centres. 
Dans les semi-conducteurs de typep, silicium ou 
germanium, la bande de valence possède un maxi- 
…. mum unique situé au centre de la zone de Brillouin. 
. Le niveau d’impureté fondamental est alors non 
dégénéré, mais la fonction d’onde ne possède pas de 
forme simple, ce qui complique grandement le pro- 
blème. Par ailleurs, il existe aussi des niveaux 
excités assez voisins de l’état fondamental. 

1 Toutes ces éventualités, présentent un caractère 
… commun : l'existence de niveaux excités voisins du 
“ niveau fondamental. Quand les centres d’impureté 
sont suffisamment rapprochés, ces niveaux donnent 
naissances à des bandes. Les fonctions de Wannier 
associées aux bandes excitées sont orthogonales 
aux fonctions de Wannier de l’état fondamental. 
Cette remarque suggère, qu’au cas où deux élec- 
trons doivent être placés sur le même ion, il puisse 
« être plus intéressant du point de vue énergétique 
nn d'attribuer au deuxième électron une orbite corres- 
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pondant à un état excité de l’ion, de façon à réduire 
l’énergie de répulsion coulombienne. 11 semble donc 
que l’on soit obligé de considérer l’ensemble des 
bandes provenant des niveaux dégénérés et quasi- 
dégénérés ; l’hamiltonien prend par conséquent 
une forme compliquée, puisqu'il fait intervenir des 
couplages entre les différentes bandes. 

- [1 semble toutefois que les principes qui nous 
ont guidé dans l’étude de modèles simples puissent 
s’appliquer à ces structures compliquées. 

Les transformations utilisées permettraient 
encore de passer continüment de la représentation 
de Bloch à celle d’Heitler et London. 

Naturellement, on aboutirait à des conclusions 
similaires : passage de l’état métallique à l’état 
isolant quand la densité électronique du système 
décroit et existence d’une énergie d’excitation e(T) 
variable avec la température. On pourrait encore 
définir des pseudo-niveaux d’énergie, variables avec 
la température. Les excitations se décomposeraient 
encore en excitations individuelles et excitations 
collectives : plasmons et spinons. 

Je remercie M. le PT Y Rocard qui a bien voulu 
accueillir dans son laboratoire et je veux expri- 
mer toute ma gratitude à M. le PT P. Aigrain qui 
m'a guidé pendant ces dernières années et dont la 
bienveillance et les conseils m'ont toujours été si 
précieux. 

Manuscrit reçu le 2 mars 1959. 
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MESURES DE COINCIDENCES RETARDÉES 
DANS LE Hô$n 


Par À. MEessen et L. GRENACS, 


Centre de Physique Nucléaire, 
Université de Louvain, Louvain (Belgique). 


Les mesures de corrélations angulaires effectuées 
dans notre laboratoire [1], [2] sur la cascade (0,410 +) 
(2,150 y) du!6Sn provenant de l’H6In ont montré une 
différence pour une source solide et une source liquide. 
Le temps de résolution du circuit de coïncidence était 
2,5.10—8 sec. Le résultat de Scharenberg [3] pour une 
source solide, avec un temps de résolution de 1.108 sec 
se rapproche nettement de notre résultat pour source 
liquide. Il était donc évident qu’une perturbation 
s’introduisait quand on admettait des gammas qui sont 
émis avec un retard supérieur à 108 sec. 

Les mesures de coïncidences retardées ont donné 
des résultats particulièrement intéressants. La cascade 
de référence (0,136 y) (0,482 y) du '8ITa, provenant 
du Z1S$Hf, indiquait un demi-temps de vie 
Tuyo = (1,0 + 0,1) 108 sec en accord avec les résul- 
tats d’autres auteurs. La droite correspondante (fig. 1) 


NT. 
Rd | 
oe se 


= —n ñ 1: 


F1G. 1. — Courbe de coïncidences retardées 
pour le 5Co, le 18H et le X6]n. 


remonte au sommet de la courbe des coïncidences 
promptes du 52Co, ce qui montre qu’il n’y a qu’une 
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seule cascade retardée qui est en Jeu. Par contre 
pour la cascade (0,410 y) (2,150 y) du 65n, on trouve 
une droite correspondante à un demi-temps de vie 
Tue = (1,0 + 0,1) 108 sec, mais qui ne remonte pas 
au sommet de la courbe des coïncidences promptes. Les 
différents points expérimentaux de. la figure corres- 
pondent à des mesures entièrement distinctes. Il y a 
deux cascades en jeu : une cascade prompte et une 
cascade retardée. Le rapport gamma retardé/gamma 
prompt est de 1,7. 

Cette mesure met en évidence qu'il faut considérer 
au lieu de la simple cascade (a) de la figure 2, l’un ou 


F1G. 2. — Schémas pour les cascades mises en évidence. 


l’autre des schémas (b) et (c) de la figure 2, avec un 
dédoublement du gamma de 2,15 MeV. Il n’est pas 
encore possible, cependant, de conclure que le second 
niveau est celui de 2,21 MeV du !16Sn provenant 
du ©6Sb. D’autres mesures seront entreprises. 

Nous remercions le D' Jones et M. Goedertier, pour 
les suggestions et discussions et pour la permission d’at- 
tirer déjà l’attention sur un de leurs résultats. Nous 
remercions aussi M. le Professeur de Hemptinne, 
Directeur du Centre, qui nous a beaucoup encouragés 
dans ce travail. 


Lettre reçue le 1er juin 1959. 
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SUR L'EFFET 
D'UN MOMENT DIPOLAIRE ÉLECTRIQUE 
DU NOYAU SUR LES NIVEAUX D'ÉNERGIE 
D'UN ATOME DANS UN CHAMP 
ÉLECTRIQUE UNIFORME 


Par JPBARRAT, 


Laboratoire de Physique de l’École Normale Supérieure, 
Paris. 


Un grand intérêt s’attache actuellement à la déter- 
mination d’une limite supérieure pour le moment dipo- 
laire électrique des particules élémentaires. On a pu 
fixer une limite très petite dans le cas du neutron 
(< e X 5.10-20 em) [1], mais les déterminations dans 
le cas du proton et de l’électron sont beaucoup moins 
précises (— e X 10-% cm) [2]. La grande finesse des 
raies de résonance magnétique observées dans les états 
fondamentaux de certains atomes (métaux alcalins, 
isotope “’Hg) [3] suggère de tenter d'établir une telle 
limite avec une grande précision en étudiant l'effet 
Stark de ces niveaux. On s’attendrait en effet à ce que 
l’interaction d’un moment dipolaire électrique nu- 
cléaire d — d1 (1 — spin nucléaire) avec un champ 


“ électrique uniforme Æ donne lieu à un effet Stark 


linéaire (interaction — d.E) que l’on essaierait de 
mettre en évidence en le séparant de l'effet Stark qua- 
dratique normal (d’ailleurs très faibles pour ces ni- 
veaux). Salpeter [2] a fait remarquer que l'effet Stark 
linéaire de l’atome d'hydrogène reste en réalité nul, 
même si l’électron porte un moment dipolaire élec- 
trique, tout au moins dans une approximation non 
relativiste. Le but de_cette note est de généraliser ce 
résultat au cas d’un atome quelconque, dans le cas où 
le noyau serait porteur d’un moment dipolaire élec- 
trique d, et de l’interpréter physiquement. 

Soient r,, ... r les coordonnées des Z électrons de 


l'atome. Le champ créé par l’électron 7 au niveau du 


noyau est 


ire AO 2 est le potentiel coulombien dû au noyau. 
1% 


2 . , 
L'hamiltonien en l’absence de champ extérieur est, 
dans une approximation nonrelativiste qui néglige le 


couplage spin-orhite. 


n, = À [Pi erv)] 


j=a L2m 


Z d Z , 

+ ZX Vlr) PS2 IV: Vlri) () 
ed et 
Pi 


où 


e? 


= le ie, hi V.. 2 
ri —7;l’ Pi 20 % ( ) 


V{ri) = 


L’interaction avec un champ électrique uniforme ZE 
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dérivant du potentiel ®(r) correspond à l’hamiltonien 
perturbateur 


Z Z 
H, = — Ÿ, e D(r,) — di.E = — Y e D(r) 
i=1 CA 
CRE 
+ Z DIV; O(r)- (6) 
ei 
puisque V; D(r;) — — E est indépendant de c. 
L’hamiltonien total s'écrit : 
H Sr H; = He + Ha (4) 
où 
/ Z Die Z 
H, = Y — eV(r,) — e Dir;) | + Y Viri) (5) 
| re AGE 
>) 
AN 
Ha Z D "ol i V(r:;) eV; D(r.,)) (6) 
\ dE : 


IT, est l’hamiltonien de l’ensemble des électrons, 
H4 celui du moment dipolaire nucléaire. h 

Le calcul se.poursuit alors comme dans le cas exa- 
miné par Salpeter [2]. On suppose que Z1, a été diago- 
nalisé et on traite 7/4 comme une perturbation ; seuls 
les termes du 17 ordre sont linéaires en d et E. On 
remarque que : 


He LS Lpi, Vlr) + ®(r 
QI Zi 2, \ -D;, r,) JE a) 
A 
: > pi, Fr) € Œir;)] () 
CR ,j=1 


puisque p; commute avec V(r,) 


+ ®(r;) pour j AZ 1. 
D’après (5), (7) s'écrit : 


HA RS re mr pl arm 
NE: ® A mnt ES en 
i>k 


Or p; commute avec p;, et avec V(r;x) sauf pour 
==) our =" #6. :Done 


—_ id BUS Go pl 
Ha = Zen À, [E:pi, He] + Zen, 2, | Di, V{ri)] 

di 
a UE 
= Ze À Pi, He] TZ, A, À Vi ij j ij 

1>j 


(9) 


Le facteur I dans le dernier terme est nul (*). Quant 
au 427 terme, ses éléments de matrice sont nuls entre 
états propres de A ayant la même valeur propre 
(c’est-à-dire ses éléments diagonaux ou entre états 


.2 
(*) D pi commute a priori avec XP & Y Viry) parce 
î j 2m it 
r] 
que D p; est lié à l’opérateur représentant une translation 


infinitésimale de l’ensemble des électrons et que 


est invariant par translation. 
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dégénérés). Donc Ha ne perturbe pas au 1er ordre les 
niveaux d'énergie de A. Le calcul est analogue et la 
conclusion reste la même si ce sont les électrons qui 
portent des moments dipolaires électriques. 

L'origine physique du résultat est l'interaction 
avec E du moment dipolaire électrique induit dans 
l’atmosphère électronique par le moment permanent 
du noyau. En effet soient 77e et ao les hamiltoniens /e 
et Ja en l’absence de champ extérieur, He, et Ha, les 
termes de 4, et de Ha dépendant du champ extérieur. 
On peut retrouver le résultat précédent en diagona- 
lisant d’abord Ho + Ha et en calculant l’effet sur ses 
états propres (déterminés au 17 ordre en d) de la per- 
turbation le, + Ha. Soient uw, les états propres 
de Hco, d'énergie Æ, (les états de même n et de « 
différents sont dégénérés et on suppose qu'ils ont la 
même parité) et One — Une + Wna les états propres 
de eo + Hao au 1% ordre en d. Ha étant impair n’a 
pas d'éléments diagonaux et : 


Wna — 


5, <Un8 a TE 7 mp (10) 
nes pe 


Dans (10), les seuls termes non nuls correspondent à 
des états u»p de parité opposée à celle de l’état ue. 
On montre comme ci-dessus (équations (7) à (9)) que : 

HUE 


es 


A PDES ILES 
So Ze ;=: 


pi, Hal 


D'où + 


LUE 
A (E Ip.) Uno 


i=1 


no — 


Doi ce S 1 ); (11 
ER MNZer: Pi ] Uno ) 


=1 


Vrax N’A pas une parité bien déterminée, puisque Wye 
et uns Sont des parités opposées. Un moment dipolaire 
est donc induit dans l’atmosphère électronique (v,, se 
déduit de u,, par une translation). En présence du 
champ E, le déplacement d’énergie linéaire en ÆE 
s'obtient en diagonalisant A, + Ha, pour les états ve 
qui peuvent rester dégénérés. Au 1€ ordre en d : 


a Dix [ar La He1| ha’ 2 = 


Z 
ee Li, I | 
EU , 


ro an Ua 


Une 


Fe die 
= < Ua [Ha;] Uno’ 22 RS < Uno [LV P(r;)| Uno’ + 


i-=1 


(12) 


Les deux termes de (12) sont égaux et opposés. Le 
1er correspond à l’effet Stark linéaire du moment dipo- 
laire d2, le 2e à l’effet Stark linéaire du moment induit 
dans l’atmosphère électronique, qui est donc égal et 
opposé au précédent. L’effet linéaire reste donc nul. 
Cette conclusion n’est valable que dans une approxi- 
mation non relativiste où l’on néglige le couplage spin- 
orbite. Des études sont en cours pour évaluer l’ordre 
de grandeur des effets linéaires qui pourraient sub- 
sister si l’on en tient compte. 


Lettre reçue le 4 juin 1959. 
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SUR UNE SOLUTION EXCEPTIONNELLE 
D'UN OSCILLATEUR A RÉLAXATION 
DU GENRE TRIODE 


Par L. SIDERIADES, 
CN ER MS MCNRAS ATEN ArdenNersele 


En électronique, la triode montée en oscillatrice 
constitue un exemple typique d’un oscillateur à rela- 
xation. L’étude de la stabilité des oscillations montre 
le rôle prépondérant de la caractéristique non li- 


néaire (@%) [1]. Si on définit le plan de phases | 


par æ —v (différence de potentiel aux bornes du cir- 
cuit bouchon dans la plaque) et y — dv/dt, on peut 
faire correspondre certaines zones de ce plan à certains 
arcs de la caractéristique. À di, = 0 il correspond 
dans le plan de phases une zone passive délimitée par 


deux droites parallèles à y’y ; et la zone active rési-. 


duelle correspond à di, 5 0. En langage topologique, 
le point singulier des courbes intégrales — qu’on prend 


pour origine — est un foyer répulsif pour la zone : 


active, attractif pour la zone passive, dans le cas des 
oscillations stables. D’après la théorie de H. Poin- 
caré [2], il existe un « cycle limite » stable qui est la 
seule courbe fermée du plan des phases et qui corres- 
pond aux oscillations effectivement observables (cf. 
Hoi). 

AU 
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N 
N Zone passive 


Zone passive 


le limite 
stable 
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Proposons-nous alors de rechercher dans quelles 
conditions il serait possible de mettre en évidence 
l'existence d’un deuxième cycle limite associé au cycle 
stable. D’après un théorème de H. Poincaré sur les 
cycles limites, ce deuxième cycle doit, être instable. 
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Considérons le schéma de la figure 2, comportant en 
plus des schémas classiques, une résistance R’ dans le 


HTC 


circuit de grille. On fait sur &,(v,) l’approximation sui- 
vante : 
5 (u 0/0 + Uo)? | Up — V5 
F 26 , P 


a, coefficient d'amplification de la triode. 

e, valeur minima de la résistance interne. 

Uo, différence de potentiel entre les potentiels de 
grille correspondant au cut-off et à la saturation. 

— #9, Valeur du potentiel de cut-off relativement à 
la masse. 

U = Ur + Wo. 

V,5, haute tension appliquée sur la plaque. 

Et on tient compte du courant de grille suivant la 
relation : 


= P'Ug + to 
po —1/p', valeur de la résistance interne grille- 
cathode. 
Cette relation n’est valable que si & est positif ; 
elle est autrement nulle. 
Dans ces conditions, tous calculs effectués, on aboutit 
au système dynamique suivant : 


dæ dy  dz d 
CR OT Z 
avec 
247 
ANT y az ] 
TEE + —— Z 
C Ê THEN 
17m : 
/, = fetet 4) 


en ayant pris soin de poser : 
Le 
r+P 
Vro, polarisation fixe de grille. 


x —1-—p'R' r! — 


B = Véo 


Le problème semble tributaire d’un espace à trois 
dimensions. Mais en fait, une projection sur le plan 


La 
my x — Ri, 


Ps tartp=0. 
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permet de retrouver la discussion plane de la stabilité, 
On peut ainsi mettre en évidence l’existence d’un deu- 
xième cycle instable dans le cas où la résistance R’ est 
grande (supérieure à la résistance interne grille- 
cathode), et où la polarisation fixe de grille V,, est 
voisine de la masse. Il se superpose alors à l’ancienne 
zone active, une zone passive plus étroite, comme 
Pindique la figure 3. Et les valeurs des paramètres « 


Fnabe 


et 8 sont telles (!) que le point singulier est attractif ; 
on peut vérifier par exemple que les conditions de sta- 
bilité de Routh-Hurwitz sont observées. On arrive 
ainsi à la conclusion suivante : 

— l’origine est un point singulier attractif dans une 
certaine zone qui l’entoure ; toute condition initiale 
sous forme d’une tension de faible amplitude conduit 
au blocage des oscillations ; | 


f 


HIc…% 


— mais il existe une zone active donnant lieu à 
l'existence d’un cycle limite stable ; on peut mettre en 
évidence l'existence de ce cycle sous la forme d’oscil- 


{) &æ = — p'R'et B = Vo — R’p'\wo — Po). 
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lations stables obtenues à partir de conditions initiales 
correspondant à une tension de grande amplitude ; 

— entre le point singulier et le cyele stable, il 
existe un cycle instable séparant les oscillations 
stables des oscillations amorties. 

En prenant une lampe dont la zone de « freinage » 
(passive), due à l’existence du courant de grille, est 
relativement faible, on peut mettre ce résultat en évi- 
dence sur l’écran d’un oscillographe cathodique repré- 
sentant le plan de phases (fig. 4). 
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RÉSULTATS DE MESURES DE TEMPS 
DE RELAXATION TRANSVERSE 7, 
PAR ÉCHOS DE SPINS 


Par Jonel Soromon, 
G. E. N. de Saclay, Gif-sur-Yvette (Seine-et-Oise). 


Nous avons mesuré le temps de relaxation trans- 
verse 7, de l’eau à 28 MHz, à une température de 24 0C 
par la méthode des Échos de Spins [1, 2], perfec- 
tionnée par Meiboom et Gill [3]. 

Pour éprouver les possibilités de l’appareillage dans 
. la mesure des temps de relaxation longs, nous avons 
mesuré les temps de relaxation d’un échantillon de ben- 
zène désoxygéné. 


TE 
115 


2 ? 


I À 


N°89, 


Ces résultats montrent que les inhomogénéités du 
champ magnétique, ainsi que la diffusion sont complè- 
tement compensées. Les défauts auxquels la méthode 
reste sensible sont les fluctuations rapides du champ 
qui, dans notre cas, obligent à faire des séquences 
d’échos très rapprochées (300 à 2 000 échos par | 
seconde). 

Les résultats des mesures pour un échantillon d’eau 
contenu dans un tube de plexiglass sont : 

Eau non dégazée : 


T2 CE RUTIRSEC 

DA? ER INSEC 
Eau désoxygénée : 

114 3,3 + 0,1 sec. 


Th = 2,65 + 0,1 sec. 


Ces valeurs ont été trouvées indépendantes du pH. 
Ces résultats ne concordent pas avec les mesures de 
Meiboom et coll. [4]. 

L’explication proposée est l'influence des parois de 
verre dans un échantillon de petite taille. En effet, si 
l’on introduit dans le tube de plexiglass des billes de | 
verre de 2 mm de diamètre, la valeur de 7’, de l’eau | 
non dégazée passe de 2,2 sec à 1,9 sec, alors que T; ne 
varie pas, aux erreurs d'expérience près. De plus cette ! 
diminution de T, dépend du pH de l’eau employée, | 
de la nature du verre et du temps pendant lequel l’eau 
reste en contaët avec la surface du verre, et l'influence 
de ces différents paramètres rend les mesures peu repro- 
ductibles. Nous n’avons fait aucune tentative d'étude 
systématique de cette interaction surface de verre-eau 
qui doit être un phénomène très complexe. M = 
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